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INFRARED ARC SPECTRUM OF GERMANIUM 
By C. C. Kiess 


ABSTRACT 


For theoretical reasons it is known that a portion of the arc spectrum of ger- 
manium must lie in the infrared. The more intense lines in this region have now 
been recorded on the specially sensitized plates of the Eastman series with the 
grating spectrographs of the National Bureau of Standards. Most of the new 
lines have been classified as combinations between the previously known *P° and 
1P° terms of the 5s configuration with new triplet and singlet terms of the 5p 
configuration. The yellow lines observed by Lunt in spark and discharge tube 
sources are found to be lines of Ge1 and, in analogy with similar lines of C1 and 
Sir, probably owe their origin to the 5s-6p transition. 


CONTENTS 


I. Introduction 
II. Experimental procedure-_-_-- 
ay: Mego atts os cece ee 


I. INTRODUCTION 


The element germanium shares with gallium and scandium the 
distinction of having its existence foretold and its properties described 
prior to its actual discovery. In 1871 Mendeléeff announced his 
periodic classification of the chemical elements and predicted three 
new members to occupy the gaps in the table as it then existed. The 
final verification of his predictions occurred in 1886, when Winkler 
discovered a new element, germanium, with the properties ascribed 
by Mendeléeff to ekasilicon. The element is rare and only small 
amounts have been prepared in a pure state. 

Investigations m9 the spectra emitted by germanium atoms in 
various states of excitation are summarized in volume 7 of Kayser’s 
Handbuch der Spectroscopie. Since the preparation of this work in 
1927, several more papers have appeared—three papers by Lang ' and 
one by Rao and Narayan,? who describe the series of Gen, Germ, 
and Getv; a paper by Rao,’? who describes the are spectrum in the 
Schumann region and gives a revised term analysis of Ger; and a 
paper by Sibaiya,* who finds that certain lines of Ger reveal no 
hyperfine structure. 

Rao’s analysis of Ge1 describes that portion of the spectrum that is 
built up on the terms *P, 1D, and 'S of the basic configuration 
4s? 4p* of the neutral atom. That part of the spectrum resulting 
from combinations between high terms and lower terms of various 


1R. J. Lang, Phys. Rev. 30, 762 (1927); Proc. Nat. Acad. Sci. 14, 32 (1928); Phys. Rev. 34, 696 (1929). 
2K. R. Rao and A. L. Narayan, Proc. Roy. Soc. (London) [A] 119, 607 (1928). 

+R. K. Rao, Proc. Roy. Soc. (London) [A] 124, 465 (1929). 

4 L. Sibaiya, Current Sci. 6, 152 (1937). 
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excited states has not yet been described. Prominent among the 
terms of the excited states are 5s*®P and 5s'P. In analogy with the 
structures that have been found for C1r® and Sir,° we may expect 
for Ger groups of lines in the infrared and in the yellow that represent 
the combinations between the higher terms of the 4s? 4p np configura- 
tions with the 5s terms. It is the purpose of this paper to describe 
the results that have been obtained in a survey of the germanium 
arc spectrum throughout the yellow, red, and infrared regions photo- 
graphically accessible. 


II. EXPERIMENTAL PROCEDURE 


The spectrograms upon which the present investigation is based were 
secured with the grating spectrographs of the National Bureau of 
Standards. For the infrared region, the grating ruled by Anderson 
with 7,500 lines per inch was employed. For the shorter regions, grat- 
ings ruled by Wood with 15,000 and 30,000 lines per inch were em- 
ployed. For each spectral region to be recorded, the appropriate sen- 
sitized plate of the Eastman series was used. 

The fig ht source was an arc between copper or silver electrodes to 
which nail amounts of GeO, were added. The arc was maintained 
at 6 or 7 amperes from 220-volt d-c mains, and operated with a fair 
degree of steadiness during the exposure periods, which ranged from 10 
or 15 minutes for the yellow region to 3 hours for the infrared. Ex- 
posures to the iron arc were juxtaposed to the germanium exposures 
to furnish the wave-length standards required in the reductions of the 
measurements. The values of the standards were selected from the 


lists adopted by the International Astronomical Union ’” or from those 
measured by Meggers.® 


III. RESULTS 


The new wave lengths of Get are presented in table 1, in which 
each entry represents the mean of two to six measurements. The 
infrared lines are sharp, only a few presenting a slightly hazy or 
widened appearance. This latter characteristic is indicated in the 
table by the letter h. On the other hand, the yellow lines are diffuse, 
many of them exhibiting an asymmetrical shading toward longer 
wave lengths. The intensities assigned to the lines are the usual 
visual estimates, and are not comparable between widely separated 
regions of the spectrum. 

5 A. Fowler and E.W.H. Selwyn, Proc. Roy. Soc. (London) [A] 118, 34 (1928); F. Paschen and G. Kruger, 
Ann. Physik [5] 7, 1 (1930). 
6 C. C. Kiess, J. Research NBS 21, 185 (1938) RP1124. 


7 Trans. Int. ‘Astron. Union 3, 77 (1928). 
’W. F. Meggers, J. Research NBS 14, 33 (1935) RP755. 





Infrared Are Spectrum of Germanium 


TABLE 1.— Wave lengths in the arc spectrum of germanium 





on | Inten- 


: Vyacctm! b bination 
sity vac Term combinatior 








11144. 8970. | 5s 3Ps—5p SP, 
11125. 8986. 5s *P3—5p *D; 
11040. 9055. 5s 1Pi—5p 88, 
11030. 9062. 5s 3P3—5p *D! 
10947. 9132. | 5s 8P7—5p 8D, 


10734. 9313. | 58 8Pg—5p 8D, 
10731. 9315. 5: 

10576. 70 9452. 

10459. 9557. | 

10405. 05 9608. 5s 3P3—5p *P, 


10382. 9628. 5s 1P?}—5p 1D, 
10305. 9701. 
10201. 9800.27 | 5s %Pi—5p*P, 
10039. 9957. | 58 $P3—5p 38, 
9957. 10039. 


9844. 10155. 
9625. 10385. 5s 3P}— 5p 5P, 
9492. 10531. 68 5s *P3—5p 'D» 
9398. 10636. 5s 8P3—5p 3P, 
9233. 10826.57 | 5p %D3—7s #P3 


9096. 5 10990. 5s 1P}—5p 183 
9068. 11023. 5s *P}—5p *P, 
8789. 11373. 5s *P}—5p 3S; 
8367. 11947. 5s *P?—5p 'Dp 
8256. 12108. 


8031. 12448, 
7833. 12762. 
7511. 13309. 5s 3Pi—5p 1S, 
6557. 15244. 5s 1P?—6p 3D, 
5967. 16752. 


5802. 17229. 5s 1P?—6p 'D» 
5801. 17233. 
5718. 17483. 
5701. 17533. 5s 8P}—6p 3D, 
5692. 17563. 5s *P}—6p *D, 


5664. 17647. 
5664. 17649. 
5656. 5 17675. 
5621. 17783. 5s *P3—6p *D, 
5616. 17801. 


5606. 17830. 


5564. 17964. 
5513. 3% : 18132. 5s *P3—6p'D; 


5265. | 18984. 
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When the wave numbers of the new infrared lines were differenced. 
it was seen that they exhibited the separations of the *P and 'P terms 
of the 5s configuration. This led to the detection of the 5p terms, as 
listed in table 2, with numerical values based on those given by Rao. 
For the present, the component *Py must be regarded as provisional] 
owing to the nonappearance of other combinations necessary to estab- 
lish it with certainty. The absence of the term 5p'P must also be 
regarded as surprising in view of the intensity with which the anal. 
ogous term appears in Sir and Ct. 


TABLE 2.—np ierms of Get 





n=5 
Term 
symbol 





Term value Term value 





3D, 17454. 25 
2D, 18723. 84 
sD, 18792. 
sP, 16832. 
sP, 17469, 
tP, 18055. 
8, 16482. 
1D, 


1p, 
1S, 


1269. 59 
69. 04 























The green and yellow lines between 5200 and 6000 A are identical 
with those previously observed by Lunt.’ In his experiments with 
sparks between metallic electrodes and with discharges through 
Geissler tubes containing GeCl,, these lines appeared under conditions 
favorable for the excitation of the neutral atom. The recent experi- 
ments at the National Bureau of Standards make it certain that they 
are emitted by neutral and not ionized atoms. Apparently they are | 
the counterparts of similar groups in Sir and C1 and represent the 
transition 5s—6p. However, the satellite lines necessary to establish | 
the multiplets with certainty are not present, and the classifications | 
given for some of them, in table 1, must be regarded as provisional. | 
The terms on which these tentative assignments are based were sug- 
gested by the occurrence among these lines of the 5s differences 250, | 
903, and 2,318 cm. 

The new terms that have resulted from this investigation belong to 
the np electron configurations. They are presented in table 2. 
According to Hund’s theory, the terms *Po, *P,, °D,, and *D, owe their 
origin to the binding of a p electron to the *Po, state of Get; the 
others, including all the singlets, have *P}, as their parent term. 


*J. Lunt, Monthly Notices Roy. Astron. Soc. 85, 38 (1924). 
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The 5p*D terms show clearly that they originate in the manner 
required by theory; the 5p*P terms exhibit this feature to a lesser 
degree. Accepting these facts as a guide, it would be possible to 
derive a complete set of terms for the 6p configuration; but for the 
present it seems more desirable to await the evidence of future obser- 
yations. 

While the work described in this paper was in progress, the wave 
lengths of some Ge 11 lines, in the yellow and red, were measured several 
times. ‘These lines appeared only at the poles of the are and were 
clearly distinguishable from the Ger lines. The mean values of these 
determinations are presented at this time, in table 3, to complete the 
observational record. 


TABLE 3.— Wave lengths of some Ge 11 lines 
Se 


4 - “4 7 - } *. . 
on “ae W “se | Term combination 








7145, 42 13991. 13 4p! 2Diy— 5p Pir 
7049. 45 | 14181. 60 4p! 2Day—5p ?Pix, 
6966. 42 14350. 62 4p’ *Diy— Sp 2Prr, 
6484. 26 15417. 71 | 5p *Pi—6s Sox, 
6336. 46 | 15777. 33 | 5p 2P34—6s 2Sox, 
6021. 09 16603. 70 | 5s "Sou —5p *Pb, 














5893. 42 5 | 16963. 38 - 58 *So%— 5p 2Phx, 








WasHineton, November 28, 1939. 
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COLORIMETRIC DETERMINATION OF ARSENIC. 
IN FERROUS AND NONFERROUS ALLOYS! 


By Clement J. Rodden 


ABSTRACT 


A method is described for the determination of arsenic and has been applied 
successfully to iron, steel, copper-base alloys, tin- and lead-base bearing metals, 
and pig lead. The arsenic is first separated from the other constituents of the 
material by distillation as the trichloride. It is then oxidized to arsenic acid and 
converted to a blue arsenic-molybdenum compound by treatment with a solution 
containing hydrazine sulfate and ammonium molybdate. The transmittancy of 
the colored solution is measured with a photoelectric colorimeter, and the amount 
of arsenic read from a curve obtained by using known amounts of arsenic. 
Alternatively, the color may be matched against standards by means of a visual 
colorimeter. ‘The method is reasonably rapid and is applicable to amounts of 
arsenic as small as 2 micrograms. 


CONTENTS 


I. Introduction______--_-__- 
II. Apparatus 
III. Weight of sample 
IV. Procedure 

V. Results 


I. INTRODUCTION 


Arsenic in ferrous and nonferrous alloys is usually separated from 
other elements by distilling it as arsenic trichloride and is then de- 
termined by (a) weighing as the sulfide, (b) converting the arsenious 
sulfide to silver arsenate and titrating the silver with potassium 
thiocyanate, or (c) neutralizing the hydrochloric acid in the distillate 
and titrating with iodine or potassium bromate.? Other methods have 
been used, especially for white metals; for example, arsenic is some- 
times separated as the element from a hydrochloric-acid solution of 
the metal or alloy by reduction with sodium hypophosphite,’ stan- 
nous chloride,‘ or calctum hypophosphite.’ The arsenic thus obtained 
is dissolved in an excess of standard iodine solution and the excess 
iodine titrated with standard sodium arsenite solution. Alternatively, 
the arsenic may be converted to arsenious acid and titrated with 
bromate.® 

! Presented —— the Division of Microchemistry at the 98th Meeting of the American Chemical Society, 
era Mass., September 11 to 15, 1939. 
BE. F. Lun ell, J. I. Hoffman, and H. A. Bright, Chemica] Analysis of Iron and Steel (John Wiley 
& ‘a Inc. .» New York, N. Y., 1931). Methods of Chemical Analyses of Metals, (American Society for 
— Materials, 1936). 
randt, Chem.-Ztg. $7, 1445, 1496 (1012). 
i. “B. Evans, Analyst ~~ 492 (1932). 


‘E. I. Fogel’son and N. V. Kalmuikova, Zavodskaya Lab. &, 584 (1936). 
‘C. W. Anderson, Ind. Eng. Chem., Anal. Ed. 9, 569 (1937). 
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For the determination of small amounts of arsenic in biological 
material, Morris and Calvery 7 recommend a modification of Denigés’ 
method,® whereby a blue complex arsenic-molybdenum compound is 
formed. In the method to be described, the formation of this blue 
color is the basis for the determination of arsenic, which has previously 
been removed from interfering elements “4 distillation as arsenic 
trichloride. Antimony, which might distill in part under certain 
conditions, does not interfere. 

The transmittancy of the colored solution is measured with 
photoelectric colorimeter (or compara- 
tor). The quantity of arsenic is read 
from a standard curve drawn from data 
obtained by treating known amounts 
of arsenic in the same manner as is de- 
scribed for the unknown. The method 
has been applied to many types of 
ferrous and nonferrous alloys. It is 
more rapid than the usual macro- 
methods and is especially useful with 
materials of very low arsenic content, 
such as pig lead. 








II. APPARATUS 


The Scherrer all-glass distilling ap- 
paratus,® shown in figure 1, was used 
to distill the arsenic trichloride. C is 
a 200-ml distilling flask, the neck of 
which is 2.56 cm in diameter. It is 
provided with an exit tube, a ther- 
mometer well, and an inlet tube for 
conveying the acids and carbon dioxide 
to the distilling flask. Both the ther- 
mometer well and the inlet tube reach 
within 3 mm of the bottom of the 
flask. The distance between the 
bottom of the flask and the exit tube leading to the condenser is 17 cm. 

The apparatus is clamped to a ring stand. C rests on an asbestos 
board (about 18 by 18 cm) having a round hole 4.5 cm in diameter. 

The photoelectric colorimeter was constructed according to the 
design of Brice * and consists of a balanced circuit using barrier-layer 
type photoelectric cells. The cylindrical absorption cells are 18 mm 
in internal diameter with the distance between the faces 2 cm. 

Examination of the colored solution with a hand spectroscope 
showed absorption in the red, and for this reason a Jena OGz filter 
which transmits from 550 my through the red was used. 

A visual colorimeter, such as the Duboscq type, may also be used 
for determining the amount of arsenic by matching against solutions 
containing known amounts of arsenic. 














Fiaure 1.—Distilling apparatus. 


71H. J. Morris and H. O. Calvery, Ind. Eng. Chem., Anal. Ed. 9, 447 (1937). 
§ G. Denigés, Pregl Festschrift, p. 27 (Mikrochemie, 1929). 

* John A. Scherrer, J. Research NBS 16, 253 (1936) R P871; 21, 95 (1938) RP1116. 
0B. A. Brice, Rev. Sci. Instr. 8, 279 (1937). Also private communications. 
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III. WEIGHT OF SAMPLE 


Nonferrous alloys, particularly bearing metals, are prone to segre- 
gate, as is shown by the variation in results of macroanalyses of 
journal bearings." The composition of the large-sized particles of 
many alloys often differs from that of the small ones, and it is usually 
inconvenient or impossible to decrease the size of all the particles to 
the point where it becomes certain that a small sample (say less than 
{0 mg) is truly representative of the material. Experiments with 
simples of brass and bronze, as prepared for macroanalyses, have 
shown that there is a sufficient difference in composition of the chips 
to cause serious error if the analysis is made on a 5- to 10-mg sample 
4 to 5 chips), but when 100-mg samples are taken, no such errors are 
encountered. 

A 100-mg sample is therefore recommended as the minimum for 
materials in the form of small chips,” unless tests have shown that 
there is no segregation, or the amount of sample is limited. For con- 
venient comparison in the cells and instrument used, the sample should 
contain at least 2 micrograms of arsenic. A 100-mg sample of com- 
mercial steel, with arsenic contents ranging from 0.005 to 0.02 per- 
cent, fulfills this condition. The same weight of sample can be used 
for cast irons. However, with high-arsenic irons (arsenic>0.03 per- 
cent) it is necessary to dilute the colored solution, as described later. 

If the arsenic content is such that 100 mg of material does not give 
sufficient arsenic, the weight of the sample may be increased. For 
example, with lead of high purity (arsenic about 0.0005 percent) a 5-g 


sample is required in order to obtain the optimum amount of arsenic. 


IV. PROCEDURE 


Carbon steel.—Place a 0.1-¢ sample in a 50-ml beaker, dissolve the 
steel in 10 ml of diluted nitric acid (1-+1), and evaporate the solution 
to dryness on a hot plate at a temperature not exceeding 130° C. Cool, 
add 5 ml of hydrochloric acid (sp gr 1.18), and transfer the solution to 
the distilling apparatus (fig. 1) with the aid of a 10-ml portion of hydro- 
chloric acid. Add 2 ml of hydrobromic acid. Mix 0.3 g of hydra- 
zine sulfate with 10 ml of hydrochloric acid (sp gr 1.18), start the 
current of carbon dioxide and transfer the mixture to the distilling 
apparatus. Insert a thermometer in the well of the distilling flask 
and provide a receiver (150-ml tall-form beaker) containing 10 ml of 
water into which the end of the condenser dips. Pass a stream of car- 
bon dioxide through the apparatus at the rate of two or three bubbles 
per second, and boil the solution gently until the temperature reaches 
il1° C. This requires approximately 15 minutes; the flask will con- 
tain about 10 ml of solution at this stage. Without interrupting the 
stream of carbon dioxide or the heating, remove the receiver and rinse 
the end of the condenser with distilled water. Add 10 ml of nitric 
acid (sp gr 1.42) to the distillate, place the beaker on a steam bath, 
and evaporate the solution to dryness. Heat the beaker and contents 
in an oven at 130° C for % to 1 hour to remove all traces of nitric acid. 
Add 10 ml of hydrazine sulfate-molybdenum solution for each 30 

"J, A. Scherrer and G. E. F. Lundell, BS J. Research 5, 891 (1930) RP237. 

'! Particle size between Nos. 18 and 40 sieves, U. 8. Standard. 

% The hydrazine sulfate-molybdenum solution is prepared as follows: Dissolve 1 g of ammonium molyb- 
date in 100 ml of 5 N sulfuric acid. Dissolve 0.15 g of hydrazine sulfate in 100 ml of water. Dilute 10 ml of 


the molybdate solution to 90 m) with water, add 1 ml of the hydrazine sulfate solution, and dilute with 
water to 100 ml. 
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micrograms of arsenic, and heat the solution on a steam bath for 10 to 
15 minutes. Cool, and dilute with hydrazine sulfate-molybdenum 
solution so that the solution contains no more than 3.0 micrograms 
of arsenic per milliliter.* This is done conveniently in 10- to 250-m| 
volumetric flasks. Measure the transmittancy of the solution with 
photoelectric colorimeter, using the hydrazine sulfate-molybdenum 
solution in the comparison cell. Read the amount of arsenic from a 
standard curve which has been drawn from data obtained by treating 
known amounts of arsenic by the same procedure as is described for 
the unknown, and correct the result for a blank carried through the 


























MICROGRAMS OF ARSENIC PER ML 
Figure 2.—Curve showing conformance to Beer’s Law. 
T =Tré@nsmittancy of solution containing arsenic. 
T.=Transmittancy of blank hydrazine-sulfate-molybdate solution. 
entire proceduré.™ The standard curve is constructed by plotting 
micrograms of arsenic with respect to transmittancy. 

The transmittancy measurements of solutions containing from 0.5 
to 3.0 micrograms of arsenic per milliliter are reproducible to about 
0.1-percent transmittancy. 

Experiments showed that the transmittancy of the solutions con- 
taining the arsenic-molybdenum compound did not change over a 
period of 24 hours. Figure 2 indicates that the colored solutions obey 
Beer’s law in the region from 0.5 to 3.0 micrograms of arsenic per 
milliliter. 

If a photoelectric colorimeter is not available, use a visual colorim- 
eter to compare the color with that obtained from known amounts 
of arsenic. A Jena OG, filter over the eyepiece of the colorimeter 
facilitates matching. 

1¢ If the approximate amount of arsenic is unknown, 10 ml of hydrazine sulfate-molybdenum solution is 
added, the solution is heated, cooled, and the transmittancy measured as directed. If the arsenic content 
exceeds 3 micrograms per milliliter, 15 ml more of the reagent is added, and the solution is again heated, 


cooled, and the transmittancy measured. This operation is repeated until the arsenic content does not 


exceed 3 micrograms per milliliter. 
18 With the reagents used, the blank ranged from 0 to 0.0002 percent of arsenic on tho basis of a 100-mg 


sample. 
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High-chromium steel.—Dissolve a 0.1-g sample in 10 ml of a mixture 
of equal parts of nitric and hydrochloric acids. Evaporate to dryness 
and proceed as with carbon steel. 

Cast iron.—Proceed as with carbon steel. 

Ferromanganese.—Proceed as with carbon steel. 

Brass, bronze, bearing metals.—Dissolve a 0.1-g sample in 5 ml of 
nitric acid (sp gr 1.42). Cool, add 8 ml of diluted sulfuric acid (1+1), 
and evaporate the solution to fumes of sulfuric acid. Cool, and add 
5 ml of water. Transfer to the distilling apparatus, and complete 
the determination as in carbon steel. 

Pig lead.—Transfer a 5-g sample to a 400-ml beaker, add 100 ml 
of difuted nitric acid (1+3), and heat on a steambath. When solu- 
tion is complete, cool, add 15 ml of sulfuric acid (sp gr 1.84), and 
evaporate to fumes of sulfuric acid. Cool, add 20 ml of water, and 
transfer to the distilling apparatus with 50 ml of hydrochloric acid 
(sp gr 1.18). Add 10 ml of hydrobromic acid. Then add 1 g of 
hydrazine sulfate to the distilling apparatus with an additional 50-ml 
portion of hydrochloric acid (sp gr 1.18). Distill as with carbon 
steel, catching the distillate in 10 ml of water contained in a 250-ml 
beaker. Distill until approximately 90 ml of distillate has passed 
over. Add 30 ml of nitric acid (sp gr 1.42) to the distillate, evaporate 
to dryness on the steambath, and complete the determination as in 
carbon steel. 


V. RESULTS 
Table 1 shows the results obtained on some National Bureau of 
Standards Standard Samples of ferrous and nonferrous alloys and on 


two samples of pig lead prepared by the American Society for Testing 
Materials for cooperative analysis. 


TaBLE 1,—Results obtained for arsenic in some ferrous and nonferrous alloys 





Arsenic Arsenic 
Materia! Material 
Present *| Found Present ® | Found 











Tin-base bearing metal, NBS 


Percent | Percent Percent ee: 
eusll “ae |4..0.H stl, NBS Standard ein 


Sample 34a ». 4 Standard Sample 54a. 


. 009 
. 009 

008 {Stainless steel (13 percent, 069 .059 }\Lead-base bearing metal, 
011 * 009 chromium), NBS Standard ; ‘061 |f NBS Standard Sample 53a. 


ee Sample 73. Bi Pret -bronze _ bearing 
012 ‘O11 
.012 { 


Ingot fron, NBS Standard < = pen = Standard 


4.0010 |\Sheet brass, NBS Standard 
Sample 37b. 


Sample 5 
- 0015 


071 O70 If Sample Ze. 4, 0018 
001 = Pig lead. 
24 |\Cast ar NBS Standard 


. 25 Sample 6d 


024 . 024 \Ferromanganese, NBS 
. 025 Standard Sample 68. 


- 068 {as iron, NBS Standard 


- 26 


4. r04 Pig lead. 























* Values reported by determinations with ne gua methods. 

» By comparing color with Duboeeq-type eolorimete 

° Analysis made on 8 6-mg sample—value obtained indie ative of nonhomogeneity of small samples. 
4 Values given to show the sensitivity and reproducibility of method. 


Wasuineton, November 1, 1939. 








U. S. DEPARTMENT OF COMMERCE NATIONAL BUREAU OF STANDARDS 
RESEARCH PAPER RP1268 


Part of Jowrnal of Research of the National Bureau of Standards, Volume 24, 
January 1940 





DETERMINATION OF URONIC ACIDS IN CELLULOSIC 
MATERIALS 


By Roy L. Whistler, Albert R. Martin, and Milton Harris ! 


ABSTRACT 


Methods for the determination of uronic acids, based on the estimation of the 
amount of carbon dioxide evolved by these acids in a given length of time during 
treatment with a boiling solution of hydrochloric acid, were found to fail when 
cellulose or other carbohydrate material was present, since these substances also 
decomposed with the evolution of some carbon dioxide. An investigation of the 
rates of evolution showed that pure uronic acids evolved carbon dioxide at a 
rapid rate and that the evolution was complete in several hours. Carbohy- 
drates, free of uronic acids, evolved carbon dioxide at a much slower rate, which 
remained approximately constant during the entire treatment with acid. It is 
shown that this difference in rates can be used as the basis for estimating uronic 
acids, even when they are mixed with large proportions of cellulosic materials. 
A detailed description of the apparatus and its manipulation is given. 


CONTENTS 


I. Introduction 
II. Description of apparatus and method 
III. Experimental results and discussion 
IV. References 


I. INTRODUCTION 


As part of a Bureau program of research relating to the structure of 
cotton fibers, it was necessary to measure the uronic acid content of 
cotton, both before and after the application of a variety of treatments. 
Such analyses necessitated the measurement of small quantities of 
uronic acids when mixed with large proportions of other carbohydrate 
material. 

Uronic acids in contact with boiling solutions of hydrochloric acid 
undergo complete decarboxylation. 


CHOH 
buon 
CHOH 











1 Research Associates at the National Bureau of Standards, representing the Textile Foundation. 
13 
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Proceeding from this observation, Lefevre and Tollens [1] * developed 
a method for the quantitative determination of uronic acids. The 
decarboxylation of the uronic acid was accomplished by boiling it in 
a solution of hydrochloric acid of specific gravity 1.06. A stream of 
carbon dioxide-free air drawn through the apparatus carried the 
evolved{carbon dioxide from the reaction flask through*a reflux con- 
denser into a Peligot tube containing water to remove traces of hydro- 
chloric acid, then through a drying tube, and finally into a weighed 
potash bulb. This procedure has since undergone numerous modifica- 
tions [2 to 11]. The subsequent methods differ for the most part in the 
manner in which the carbon dioxide is collected and determined, and 
in the length of time necessary to bring the reaction to completion. 
The length of time recommended as necessary for complete de- 
carboxylation has ranged from 3% to 10 hours. 

Although these methods gave good results when applied to pure 
uronic acids, considerable difficulty was encountered when they were 
used on such materials as cotton fibers or other carbohydrate mix- 
tures. Under the drastic conditions employed in the analysis, it was 
found that uronic acid-free carbohydrates decomposed with the 
evolution of carbon dioxide. For example, when pure glucose was 
treated with a boiling 12-percent solution of hydrochloric acid, carbon 
dioxide was continuously evolved; the evolution was appreciable even 
after 30 hours. A study of the rates of evolution of carbon dioxide 
revealed that in the case of pure uronic acids there occurred an initial 
rapid evolution of carbon dioxide which ceased after several hours. 
On the other hand, carbohydrates which did not contain uronic acid 
groups evolved carbon dioxide at a much slower but approximately 
constant rate during the entire treatment with acid. A mixture of 
uronic acids with other carbohydrate material gave an initial rapid 
evolution of carbon dioxide followed by a slower but steady evolution 
of the gas. 

During the time that the present work was in progress, there ap- 
peared three papers [12, 13, 14] from other laboratories, reporting 
findings relative to the evolution of carbon dioxide from a number of 
pure sugars and polysaccharides. Norman [13] pointed out that the 
presence of uronic acid groups could be inferred from the rate of 
evolution of carbon dioxide (a fact which we had noted independently) 
and stated that no quantitative significance could be attached to the 
small yields of carbon dioxide from polysaccharides. 

In the present investigation, a critical study was made of the 
factors which influence the rates of evolution of carbon dioxide from 
uronic acids, glucose, and cellulose fibers during treatment with 
hydrochloric acid. As a result of these studies, it has been possible 
to develop a new method for the precise estimation of small quanti- 
ties of uronic acids in the presence of cellulosic materials. 


II. DESCRIPTION OF APPARATUS AND METHOD 


The apparatus which was developed in this investigation is shown 
in figure 1. Two assemblies of the type shown were employed 
simultaneously; the same source of carrier gas and the same oil bath 
were used for both. 


2 Figures in brackets indicate the literature references at the end of this paper. 
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Nitrogen, which is used as the carrier gas for the evolved carbon 
dioxide, enters the apparatus through an empty safety bottle, A. 
It next passes through an alkaline solution of pyrogallol, B. The 
inlet tube in this bottle is drawn out to a small orifice which produces 
yery fine bubbles. From B the gas passes through two absorption 
towers, C, filled with soda lime, into a second safety bottle, D, which 
is provided with a mercury manometer, E, and then enters the 
reaction flask, F', by way of a side arm whose outlet is 10 to 15 mm 
above the surface of the liquid in the flask. The size of the flask 
depends upon the type of material and the size of the sample to be 
analyzed. In most of the experiments a 500-ml reaction flask was 
used. From the reaction flask the gas passes through a 40-cm 
reflux condenser, G, and into a bubbling tower, H, containing 
approximately 60 ml of concentrated sulfuric acid. The sulfuric 
acid serves to remove interfering decomposition products which are 
carried over from the reaction flask. The gas next passes through 
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FigurRE 1.—Apparatus for determination of the rate of evolution of carbon dioxide 
from uronic acids or materials containing uronic acids during treatment with 
hydrochloric acid. 


(The apparatus is not drawn to scale. See text for dimensions.) 


the U-tube, J, which is filled with anhydrous copper sulfate,* then 
through the tube, J, which contains phosphorus pentoxide, and 
finally through the carbon dioxide-absorption tube, K, containing 
ascarite, backed by phosphorus pentoxide. The absorption tubes, 
K, are connected into the train by means of mercury-cup seals [15]. 
This type of connection makes possible a rapid exchange of the ab- 
sorption tubes. Tube K is protected by a soda-lime tube, L, which 
is followed by a calibrated flowmeter, M, for estimating the rate of 
flow of nitrogen through the apparatus. 

The reaction flask is immersed in a vessel containing about 16 
liters of hydrogenated cottonseed oil. A bath temperature of 130° C 
was Pes to be optimum for maintaining a steady but gentle boiling 
of the reaction mixture. The bath is brought to the operating 
temperature by means of two electric immersion heaters, one of 500 
and one of 1,000 watts. When the desired temperature is reached, 

} An aniline-hydrochloric acid mixture may also be used to remove volatile products, but in our work 


it was found to be less efficient than H2S04. é 
‘A precaution necessary in systems where chlorine or hydrogen sulfide might be formed. 
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the 500-watt heater alone is sufficient to maintain thermal constancy 
within +0.2° C. The time required to raise the temperature of the 
bath to 130° C is approximately 50 minutes. 

The sample to be analyzed is placed in the reaction flask, F. The 
optimum size of the sample depends upon its uronic acid content, 
When working with a material such as cotton, which contains g 
relatively small proportion of uronic acids, it is advisable to use samples 
of approximately 10 g, although good results may be obtained with 
smaller samples. A sample of this size is readily covered by 300 ml 
of the hydrochloric acid solution. Varying the ratio of acid solution 
to sample within reasonable limits does not affect the results. Thorough 
wetting of the sample by the acid solution is accomplished by evacuat- 
ing the air in the flask with a filter pump. 

The flask is placed in position in the oil bath so that the oil level is 
3 to 4 mm lower than the liquid level inside of the flask. This pre- 
caution is taken to prevent the baking of small bits of sample which 
may be splashed against the sides of the flask. Nitrogen, at the rate 
of about 10 liters per hour, is passed through the apparatus until 
the ascarite tube, BR shows no further gain in weight. This operation 
requires about 30 minutes, during which time the temperature of the 
oil bath is slowly raised to 50° C. When the apparatus is free of 
carbon dioxide, both heating units are turned on and the temperature 
brought to 130°C. This procedure is always carefully followed in 
order to assure the same preliminary heating for all samples. The 
point of zero time is taken as the time at which the bath reaches 
130°C. At that time the ascarite tube, K, is removed for weighing 
and a second weighed ascarite tube inserted in its place. At the end 
of 1 hour the second tube is removed for weighing and replaced by the 
first. This process is repeated at intervals of 1 hour for the duration 
of the analysis. When analyses are made on pure uronic acids or 
materials rich in uronic acids, a small amount of carbon dioxide is 
evolved by the time the temperature of the bath reaches 130° C. 
In these cases, this amount is measured and added to that evolved 
during the first hour. 

Previous investigators have suggested that the tube entering the 
reaction flask, F, extend into the reaction mixture in order to prevent 
bumping. However, it was found that there is no tendency to bump 
when a 12-percent solution of acid and an oil-bath temperature of 
130° C are used. On the other hand, it is found advantageous to have 
this tube terminate above the surface of the liquid. With this arrange- 
ment there is no tendency toward different amounts of agitation re- 
sulting from slight variations in the rate of flow of the carrier gas. 
Furthermore, there is less tendency for the gas stream to carry inter- 
fering products from the reaction mixture into the sulfuric acid tower, 
H. That no secondary reaction occurs from the action of sulfuric 
acid on the volatile substances, which interferes with the determina- 
tion, was demonstrated by substituting a U-tube cooled in a mixture 
of solid carbon dioxide and trichloroethylene for the sulfuric acid 
tower, H. Identical results were obtained under both conditions. 

The length of time required to sweep the carbon dioxide from the 
reaction flask into the absorption tubes was determined by the use 
of a known quantity of pure sodium carbonate. The sample was 
weighed in a glass tube, and the ends were sealed with paraffin. The 
tube was placed in an acid solution in the reaction flask and the train 
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swept out with nitrogen. The temperature of the bath was then 
raised until the paraffin plugs loosened and allowed the acid to come in 
contact with the sample of carbonate. Since the reaction is extremely 
rapid, it was easy to determine the moment at which the carbon 
dioxide was released into the system. This procedure was repeated 
several times, a different rate of nitrogen flow being used in each case. 
It was found that with a gas flow of 10 liters per hour, all the carbon 
dioxide was absorbed by the ascarite in less than 30 minutes. 


III. EXPERIMENTAL RESULTS AND DISCUSSION 


The rate of evolution of carbon dioxide from pure galacturonic acid 
monohydrate in 11-, 12-, and 13-percent solutions ° of hydrochloric 
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FiguRE 2.—Rate of evolution of carbon dioxide from galacturonic acid monohydrate. 





acid is shown in figure 2. The curves show that an initial, rapid evolu- 
tion of carbon dioxide occurs and that the evolution ceases after about 
‘hours. <A total of 208 mg of carbon dioxide per gram of galacturonic 
aid monohydrate was obtained (the theoretical yield is 207.5 mg). 
Increasing the concentration of acid increases the rate of evolution 
but does not affect the total amount of gas evolved. 

Similar treatment of pure glucose in solutions of hydrochloric acid 
tives very different results, as shown in figure 3. The amounts of 
carbon dioxide evolved are small * compared with the amounts obtained 
rom galacturonic acid under similar conditions. The rate of evolu- 
tion increases with increasing acid concentration; but unlike the case 
of the galacturonic acid, the evolution of gas continues, except in the 


' These solutions contained 110, 120. and 130 g of hydrogen chloride per liter of solution. 
‘ Because of the small amount of CO: evolved, 10-g samples of glucose were used for each analysis. 
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most concentrated acid solutions, at an approximately constant rate 
for the duration of the analysis. 
The amount of carbon dioxide obtained from pure glucose is very 
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Fiaure 3.—Rate of evolution of carbon dioxide from glucose. 


small compared with that obtained from pure galacturonic acid; but 
it is obvious that in naturally occurring cellulosic materials, such as 
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TIME IN HOURS 
Ficure 4.—Rate of evolution of carbon dioxide from dewaxed cotton. 





cotton, which contain only small amounts of uronic acids, the amount 
of carbon dioxide evolved by that portion of the material which does 
not contain uronic acid groups may be an appreciable part of the total. 
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This is shown by the curves for dewaxed cotton’ in figure 4. An 
examination of these curves indicates that they are probably compos- 
ites of curves of the types found in figures 2 and 3. The rapid evolu- 
tion of carbon dioxide during the first 3 hours, as indicated by the 
steep slopes of the curves, is characteristic of uronic acids. After the 
3-hour period the evolution of gas does not cease but continues at a 
slower but constant rate characteristic of carbohydrates which do not 
contain uronic acid groups. As in the previous experiments, the rate 
of evolution of gas increases with increasing acid concentrations. 

On the basis of the results shown in figures 2, 3, and 4, a 12-percent 
solution (3.290 N) of hydrochloric acid was chosen as a satisfactory 
concentration for the decarboxylation of uronic acids in the presence 
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FiaurE 5.—Rate of evolution of carbon dioxide from glucose, galacturonic acid, and 
a mixture of both. 


of other carbohydrate material.’ In lower concentrations of acid the 
decarboxylation of the uronic acids proceeds too slowly for convenience, 
whereas in acids of higher concentration the decomposition of the 
carbohydrate itself is too rapid and erratic. Since the rate of evolution 
of carbon dioxide is appreciably affected by variations in acid strength, 
it is essential that the same concentration (within 0.02 percent) of acid 
be used in all comparative studies. 

That the results obtained with dewaxed cotton probably represent 
the behavior of a mixture of uronic acids and other carbohydrate 
materials was further indicated by studying the evolution of carbon 
dioxide from glucose, galacturonic acid, and a known mixture of the 
two. The results are shown in figure 5. Curve C, which represents 

7 Dewaxed cotton refers to cotton which has undergone 24-hour extractions each with alcohol and with 


diethyl ether. 
’ This concentration was used in all subsequent experiments. 
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the total amount of carbon dioxide evolved from a mixture of 10.00 g 
of glucose and 0.100 g of galacturonic acid monohydrate, is similar in 
shape to the curves in figure 4. The total amount of carbon dioxide 
evolved by the mixture at any time, as shown by curve C, is equal to 
the sum of the amounts evolved by each of the components in the 
same time, as shown by curves A and B. The good agreement of 
the data is seen in table 1. 

Conversely, it follows that the amount of carbon dioxide evolved 
from uronic acid groups in this or similar mixtures can be readily 
estimated from the differences between the ordinates of the parallel 
portions of the curves B and C. The average of these differences, 
taken from the data in table 1, is 20.5 mg of carbon dioxide (theory, 
20.75 mg). 


TABLE 1.—Evolution of carbon dioxide from galacturonic acid, glucose, and a 
mixture of both 
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The curves shown in figure 6 were obtained from samples of dewaxed 
cotton and from cotton which had been purified by boiling for 7% 
hours in a 1-percent solution of sodium hydroxide according to the 
methods recommended for the preparation of standard cotton cel- 
lulose [16, 17, 18]. With the exception of the induction period during 
the first 2 hours, the curve for purified cotton is similar in shape to 
that for glucose. It shows no initial, rapid rise, characteristic of the 
uronic acids, and indicates that cotton purified in this manner is free 
from a measurable amount of uronic acids. The amount of carbon 
dioxide evolved from the uronic groups in the dewaxed cotton is ob- 
tained by taking the average distance between the parallel portions 
of the curves, which gives a value of 2.4 mg of carbon dioxide per 
gram of cotton. 

It has generally been assumed that the uronic acids in native cotton 
are contained in the pectic ® substance of the fiber. Further evidence 
for this assumption is shown in figure 7. The carbon dioxide-rate 
curve of the pectin ® is very similar in shape to that obtained for 
galacturonic acid (fig. 2), and reaches a maximum at 208 mg of car- 

* Pectic substance has been defined by the ACS committee on nomenclature of _— [19] as a complex 


carbohydrate derivative which is characterized by the presence of galacturonic aci 


% The pectin used in these opetnes was obtained from the cotton by the well-known method of 
extraction with ammonium oxalate. A description of the preparation and properties of this pectin is being 
prepared for publication. 
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FigurE 6.—Rate of evolution of carbon dioxide from samples of dewaxed and purified 
cotton. 
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FicurE 7.—Rate of evolution of carbon dioxide from a sample of pectin from colton. 
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bon dioxide per gram of pectin. From this value a factor of 4.8" 
for converting carbon dioxide, formed from uronic acids, to the pectin 
extracted from cotton is obtained. The pectin content of the cotton 
used in the present work may then be roughly estimated by multiply- 
ing 2.4 mg (from the curves in fig. 6) by 4.8, which gives a value of 
11.5 mg of pectin per gram, or 1.15 percent. As is to be expected, 
this value is higher than previously reported values obtained by 
isolation methods [20, 21]. 

The uronic acid contents of a number of other cellulosic fibers were 
also determined. Table 2 shows the results obtained from decorti- 
cated ramie, retted flax, and purified samples of both. Purification 
was accomplished by treatment of the fibers with a boiling 1-percent 
solution of sodium hydroxide for 4 hours by the methods described 
for cellulose. The data for the decorticated ramie and retted flax 
are similar to those obtained for dewaxed cotton, inasmuch as they 
show an initial, rapid evolution of carbon dioxide, which becomes 
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Figure 8.—Rate of evolution of carbon dioxide from samples of rayon, glucose, 
and purified cotton 


slower but approximately constant after 4 hours. During the first 4 
hours considerably more carbon dioxide was evolved by both ramie 
and flax than was evolved in the same period by cotton. The value 
for the amount of carbon dioxide evolved by the uronic acid groups, 
obtained by subtracting the value for the purified from that for the 
unpurified fiber, reaches a maximum after about 4 hours and then 
remains constant. An average of 10.4 and 10.2 mg of carbon dioxide 
per gram of ramie and flax, respectively, was obtained. 

The results obtained from three different samples of rayon are shown 
in figure 8. The curves for viscose and cuprammonium rayon are 
nearly identical with that for glucose. A sample of acetate rayon 
evolved carbon dioxide at a somewhat lower rate. However, if the 
data for the acetate rayon are calculated on the basis of milligrams of 
carbon dioxide per gram of cellulose instead of cellulose acetate, the 

ul This factor is valid only if it is assumed that the pectin has not been appreciably altered during removal 


from the fiber. Justification for the use of such a factor will depend on obtaining similar results on samples 
of pectin removed from cotton by several different methods. Work along this line is in progres?. 
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results are very similar to those for the viscose or cuprammonium 
rayon. The curves for these samples indicate that they are free of 
uronic acids. 


TaBLE 2.—Amount of carbon dioxide evolved from the uronic acid groups of ramie 
and flax fibers 
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A carbon dioxide-rate curve for purified cotton cellulose is included 
for comparison. It is of interest to note that the slope of this curve is 
quite different from that of glucose or the regenerated fibers. 
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THERMAL EXPANSION OF SOME CHROMIUM- 
VANADIUM STEELS 


By Peter Hidnert 


ABSTRACT 


Data on the linear thermal expansion of several SAE chromium-vanadium steels 
have been published previously. The present paper gives the results of an in- 
vestigation at this Bureau on the linear thermal expansion of two other SAE 
chromium-vanadium steels (SAE 1625 and 6150) in the annealed condition and in 
the quenched and tempered condition. These steels were investigated at various 
temperatures between 20° and 650°C. The coefficients of expansion of the 
hardened sample tempered at 400°F (204°C) or 800°F (427°C) are less than the 
coefficients of expansion of the annealed sample. Tempering at 1,200°F (650°C) 
ncreased the expansion to that of the annealed sample. For the range from 20° 
to 650°C, the average coefficients of expansion of all of the samples investigated 
vary from 13.8 10-° to 14.8 10-° per degree centigrade. 


CONTENTS 


I. Introduction 
II. Materials investigated 
III. Apparatus 


V. References 
I. INTRODUCTION 


Coefficients of linear thermal expansion of SAE 6115,! 6135, and 
6140 } chromium-vanadium steels, have been determined by Souder 
and Hidnert [1, 2] ? and republished in 1939 in the Metals Handbook 
of the American Society for Metals. The present paper gives the 
results of an investigation of the linear thermal expansion of other 
SAE chromium-vanadium steels which are strictly in accordance with 
SAE specifications. 

Chromium-vanadium steels are used for a variety of purposes. 
Publications by Vanadium Corporation of America [3] and Priestley 
[4] give information about the mechanical properties and uses of these 


steels. 
II. MATERIALS INVESTIGATED 


The chemical composition and the heat treatment of five samples 
that were investigated are given in table 1. Two of these samples are 
SAE 6125 chromium-vanadium steels and the other three samples 
SAE 6150 chromium-vanadium steels. 


1 Chemical composition not strictly in accordance with SAE specification. 
?The noumbers in brackets indicate the references at the end of this paper. 
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III. APPARATUS 


The samples were investigated at various a between 
10° and 650° C with the precision comparator type of thermal-expan- 
jon apparatus damaiool by Souder and Hidnert [5]. The white 
furnace Shown at the left of figure 1 of their publication was used. 
figure 5 of the same publication shows the details of mounting the 

ample in the furnace. The distance between the two sharp V 
notches on each sample was 300 mm. 

Platinum-osmium observation wires (6% percent osmium) were 
ised. The diameter of each observation wire was 0.050 mm. A vane 
weighing 18 g was attached at the bottom of each observation wire 
and suspended in oil to damp out vibrations. 


IV. RESULTS 


The observations obtained on the linear thermal expansion of the 
ive samples of chromium-vanadium steels at various temperatures 
between 20° and 650° C are shown in figures 1 and 2. Each expan- 
sion curve was plotted from a different origin. Two expansion tests 
were made on samples 1611 and 1613, for in the first tests the observa- 
tions on cooling lie appreciably below the expansion curves obtained 
on heating. Results obtained in a second test represent data on 
material after the heat treatment incident to the first test. 

In each test except the first tests on samples 1611 and 1613, the 
observations on cooling are very close to the expansion curve. This 
occurs in each case when the sample had been annealed or if the 
sample had been heated to 1,200° E (about 650° C) after quenching. 

The average coefficients of expansion in table 1, were obtained 
from the expansion curves of figures 1 and 2. This table also shows 
the difference in length before and after each expansion test. 

A comparison of the coefficients of expansion of samples 1610 and 
1611 shows that the expansion of the hardened sample which had 
been tempered at 400° F (204° C) is less than the expansion of the 
annealed sample. The tempering produced by the first test to 1,200° F. 
(650° C) increased the expansion (in the second test) to that of the 
annealed sample. This effect is also indicated by the contraction 
curve of the first test and the expansion and contraction curves of the 
second test. 

Similarly, for samples 1612 and 1613, it is evident that the expan- 
sion of the hardened sample which had been tempered at 800° F 
(427° C) is less than the annealed sample. The tempering produced 
* heating to 1,200° F (650° C) in the first test, also increased the 
expansion (in the second test) to that of the annealed sample. This 
effect is also indicated by the contraction curve of the first test and 
the expansion and contraction curves of the second test. 

The tempering at 1,200° F (650° C) for sample 1614 is comparable 
to the heat treatment incident to the first tests on samples 1611 and 
1613 and produced the same expansion as that of the corresponding 
annealed sample. 
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Figure 1.—Linear thermal expansion of two samples of SAE 6125 chromium- 
vanadium steel. 


C, 0.22; Mn, 0.75; Si, 0.27; P, 0.019; S, 0.033; Cr, 0.96; V, 0.17 percent. 


Sample 1610—Annealed at 1,650° F for 3 hours and furnace cooled; sample 161 eae for % hour at 1,475° F, 
quenched in oi] and tempered 4 hours at 400° 





Hie 


Fi 


Hidnert] Expansion of Some Chromium-Vanadium Steels 











ke 
z 
uJ 
O 
a 
ud 
a. 
z 


- 
| 
| 

N 
° 

Z| 

o 

O 

2 

4 

a 

x< 

fi 


a 
HEATING COOLING 


FIRST TEST ° ° 


SECOND TEST 4 a 
—— ' —_ 











100 200 300 400 
TEMPERATURE 


FiaureE 2.—Linear thermal expansion of three samples of SAE 6150 chromium 
vanadium steel. 


C, 0.53; Mn, 0.80; Si, 0.15; P, 0.615; S, 0.020; Cr, 1.02; V, 0.17 percent. 
Sample 1612—Annealed at 1,575 ° F for 3 hours and furnace cooled; sample 1613—heated for % hour at 1,575° F, 
quenched in oil, heated for 3 hours at 800° F and air cooled; sample 1614—heated for 4 hour at 1,575° F, 
quenched in oiJ, heated for 3 hours at 1,200° F and air cooled. 
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SENSITIVE ANEROID DIAPHRAGM CAPSULE WITH NO 
DEFLECTION ABOVE A SELECTED PRESSURE 


By W. G. Brombacher, V. H. Goerke, and F. Cordero 


ABSTRACT 


A sensitive diaphragm capsule has been designed in which the two corrugated 
diaphragms nest into each other at external air pressures above a desired value. 
Evacuated capsules of this type have particular application in measuring pressure 
with radio sondes of the Diamond-Hinman-Dunmore type, where the deflection 
of the contact arm is fixed. At the ground level the diaphragms are designed to 
nest into each other and deflection does not begin until the air pressure is reduced 
to the value called the cut-off pressure. Several elements of two capsules each 
were constructed with a cut-off pressure of 140 millibars. When installed in a 
radio sonde in place of the usual type, which is responsive over the entire range of 
pressures, a sevenfold increase in sensitivity in pressure measurement was obtained 
at altitudes above 46,000 feet. 


In measuring air pressure with radio sondes of the Diamond- 
Hinman-Dunmore type,’ where the motion of the pressure contactor 
is a fixed distance, it 1s desirable, in some cases, to increase the sensi- 


tivity at high altitudes. To secure this increase in sensitivity, a 
disphragm capsule has been designed at the National Bureau of Stand- 
ards, the two diaphragms of which, when evacuated and sealed off in 
accordance with the usual practice for capsules used to measure atmos- 
pheric pressure, nest into each other at air pressures greater than a 
selected value. At pressures less than this cut-off value, the capsule 
deflects normally. The nesting of the diaphragms prevents damage 
to the relatively sensitive diaphragms, which might be caused by exces- 
sive differential pressures. The new diaphragm element is designed 
so that the contact arm of the radio sonde deflects over the same dis- 
tance for the pressure interval from the cut-off to zero pressures, as 
does the ordinary element from sea level to zero pressures. 

The capsule, at zero pressure, is shown diagrammatically in figure 1. 
The connections needed for evacuating the capsule and for attach- 
ment to the radio sonde are not shown in figure 1. It is apparent 
from the form of the capsule that, as the air is evacuated from its 
interior, it will gradually collapse and finally, at a definite differential 
pressure, the two diaphragms will nest into each other and remain 
nested for larger external pressures. The evacuated capsule starts 
deflecting when the air pressure is reduced to a pressure equal to the 
differential pressure at which nesting took place during evacuation. 
This is known as the cut-off pressure. 

A pressure element of this type, consisting of two diaphragm cap- 
sules installed in a Diamond-Hinman-Dunmore radio sonde made by 
‘Harry Diamond, Wilbur S. Hinman, Jr., and Francis W. Dunmore, A method for the investigation of upper- 
ur phenomenc and its application to radio meteorography, J. Research NBS 20, 369 (1938) RP 1082. 
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Friez, is shown in figure 2. Two such elements were constructed. 
The diaphragms are 0.003 inch thick and 2.5 inches in diameter. 
They were formed hydraulically from soft beryllium copper, which 
was later hardened by heat treatment. Separate dies were used 
to form the upper and lower diaphragms. The complete pressure ele- 
ments have a deflection of about 0.22 inch at the center for the pres- 
sure range from 140 to 0 millibars (105 to 0 mm of mercury). At 
greater air pressures, the diaphragms of each capsule nest into each 
other. As installed in the radio sonde, the first contact is made at 
about 140 millibars (105 mm of mercury), equivalent to an altitude of 
about 46,000 feet. Since the pressure element ordinarily installed in 
the radio sonde deflects over the range of contacts in the pressure 
range from 0 to 1,040 millibars, the gain in sensitivity on a pressure 
basis is about 7% to 1. The effect of temperature change (0° to +30° 
C) on the deflections was found to be somewhat less than that to be 
expected on the basis of the temperature coefficient of the modulus of 
elasticity of beryllium copper. 

In cooperation with H. Diamond and W. S. Hinman, one of the 
radio sondes was sent up attached to six sounding balloons, which 
reached an altitude of 85,500 feet. Seventy of the eighty contacts 
of the range operated during the flight. A temperature record was 
also secured. The minimum pressure obtained was 24 millibars (18 


EET 


Figure 1.—Capsule in which the diaphragms nest into each other as the internal 
pressure is decreased or the external pressure increased. 


mm of mercury). The action was satisfactory in all respects, includ- 
ing smooth motion of the contactor over the contacts. 

Other designs may be indicated, using the formulas developed by 
W. A. Wildhack and V. H. Goerke* Thus two beryllium-copper 
capsules similar to those described above, and with the same shape 
of corrugations but of material 0.0045 inch thick, would have the 
same deflection in the pressure range 0 to 256 millibars (0 to 200 
millimeters of mercury). 

The effect of friction at the sliding contact should be considered. 
If the friction at the point of contact and the deflection range of the 
contact arm remain constant, while the stiffness of the pressure 
element and the multiplication of the mechanism are varied, the 
differential pressure necessary to balance the friction varies directly 
as the stiffness and inversely as the square of the diameter of the 
diaphragm capsule. 

If data are required at air pressures above the cut-off point of the 
elements described, the fovol type of pressure element can also be 
used, but should be arranged so as to record only in this pressure 
range. The sensitive element contactor would be so arranged that 
it would traverse the samefcontacts at pressures below its cut-off 
pressure. 

This diaphragm capsule was developed at the National Bureau of 
Standards as a part of the program of upper-air investigations under- 
taken in cooperation with the National Geographic Society. 


Wasuineton, November 27, 1939. 


2 Corrugated Metal Diaphragms for Aircraft Pressure Measuring Instruments. Nat. Adv. Comm. 
Aero. Technical Note 738 (1939). 
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Figure 2.—New type pressure element installed in Diamond-Hinman-Dunmore 
radio sonde. 


A is the contact arm. 
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PHYSICAL PROPERTIES OF SOME PURIFIED ALIPHATIC 
HYDROCARBONS 


By Donald B. Brooks, Frank L. Howard, and Hugh C. Crafton, Jr. 


ABSTRACT 


In an investigation of the suitability of various paraffin hydrocarbons as 
constituents of aviation fuel, which is being conducted for the National Advisory 
Committee for Aeronautics, the Navy Bureau of Aeronautics, and the Army 
Air Corps, four olefin and seven paraffin hydrocarbons have been obtained in a 
state of high purity. Eight of these materials were synthesized, one was isolated 
from & commercial synthetic crude, and two were obtained from commercial 
sources. All were purified by distillation in automatically controlled fraction- 
ating columns of high efficiency. 

The measured physical properties of these’ materials included freezing point, 
boiling point and its variation with pressure, refractive index and density and 
their variations with temperature. 


CONTENTS 


I. Introduction 
II. Fractionating columns 
III. Preparation of materials 
1. Preparation of 3,3-dimethylbutene-1, 2,3-dimethylbutene-1, and 
2,3-dimethylbutene-2 
. Preparation of 2,2-dimethylbutane__-_.....__..__-_..-._____- 
. Preparation of 2,3-dimethylbutane 
. Preparation of 2,3,3-trimetbylbutene-1 and 2,2,3-trimethy!- 
butane 
. Preparation of 3-ethylpentane 
. Isolation of 2,3,4-trimethylpentane 
. Purification of n-heptane and 2,2,4-trimethy!pentane________. 
IV. Determination of physical properties 
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I. INTRODUCTION 


The materials described herein were prepared in connection with 
an investigation of the suitability of various paraffin hydrocarbons as 
constituents of aviation fuel. This investigation is being conducted 
at the National Bureau of Standards for the National Advisory Com- 
mittee for Aeronautics, the Navy Bureau of Aeronautics, and the 
Army Air Corps. In the previous paper [1]! of this series, the 
preparation and properties of purified 2,2,3-trimethylpentane were 
described. 

Eight of the hydrocarbons described herein were synthesized, one 
was isolated by fractionation from a commercial synthetic crude, and 
two were obtained from commercial sources. All of these materials 
were further purified by distillation in automatically controlled frac- 
tionating columns of high efficiency. 


' Figures in brackets indicate the literature references at the end of this paper. 33 
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II. FRACTIONATING COLUMNS 


The high purity of the compounds described in this paper resulted 
primarily from the efficiency of the two fractionating columns used 
for the final purification. Both columns operate automatically and 
er throughout a distillation, whether it takes 4 days or 4 
weeks. 

Both columns are built entirely of Pyrex glass. The smaller is g 
150-plate bubble-cap column, the fractionating elements of which 
follow the design of Bruun [2]. The other column is a 600- by 4-cm 
helix-packed column. Both columns have silvered vacuum jackets, 
heating elements wound around the latter, and 1 inch of asbestos 
insulation outside the heating elements. Five-liter still-pots are sealed 
on the columns. The head of each column has a thermometer well. 
The columns operate under total reflux with intermittent sampling 
automatically timed. 

Other columns similar to the bubble-cap have shown efficiencies of 
over 80 theoretical plates for a 120-plate column; therefore, it is be- 
lieved the efficiency of the 150-plate bubble-cap column is close to 100 
theoretical plates. By comparison of the distillation of the same or 
similar materials in the two columns, it was evident that the 600- by 
4-cm packed column was more efficient than the bubble-cap column. 

The criteria which guided the selection of samples for measurement 
of physical properties were the constancy of boiling point and the 
refractive index of the fractions, and other factors of the particular 
run which experience had shown might affect the purity of the fraction. 


III. PREPARATION OF MATERIALS 


1. PREPARATION OF 3,3-DIMETHYLBUTENE-1, 2,3-DIMETHYLBU.- 
TENE-1, AND 2,3-DIMETHYLBUTENE-2 


Flow diagram: 
HeCl, 


Acetone (Ref 3) 





(Ref 3) 





> Mg amalgam 


Mg 
. (Ref. 4) 
pinacol hydrate —> 





. Na ; dehydrate 
pinacolone > pinacolyl alcohol > 
distill 
hexenes —~ 











3,3-dimethylbutene-1 (51%) 
2,3-dimethylbutene-1 (40%) 
2,3-dimethylbutene-2 (9%) 


Al: (SO,)s 
: catalyst 
Excess 3,3-dimethylbutene-1 
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3,3-dimethylbutene-1 (7%) 
2,3-dimethylbutene-1 (28%) 
2,3-dimethylbutene-2 (65%) 


A. Pinacol hydrate-——Pinacol hydrate was prepared from acetone 
and magnesium amalgam by a modification of the procedure given 
in Organic Syntheses [3]. 

A 5-liter three-necked flask was fitted with a 500-ml dropping 
funnel and a 150- by 2-cm condenser in the two smaller necks. Water 
at 5° to 10°C was circulated through the condenser jacket. The 
large center neck was stoppered with a cork. In the flask were placed 
120 g of magnesium turnings and 500 ml of benzene. 

(Technical benzene and acetone for this preparation were dried for 2 
days over calcium chloride and decanted. No increase in yield was 
found when these materials were further purified.) A solution of 135 ¢ 
of mercuric chloride in 260 ml of acetone was added all at once, and 
the very vigorous reaction controlled by adding portions of 700 ml of 
benzene as needed. In this way all the magnesium was amalgamated 
equally before the magnesium pinacolate had swelled up to the point 
where addition of mercuric chloride caused local amalgamation and 
separation of free mercury. After the benzene had been added, 500 
ml of acetone was added slowly, and then a solution of 300 ml of 
benzene and 390 ml of acetone was added. The swollen magnesium 
pinacolate was worked down to the bottom of the flask, when neces- 
sary, by inserting a wooden paddle in the large neck of the flask and 
stirrmg. The product was worked up in the way described in Organic 
Syntheses [3]. Concentration of the pinacol hydrate solution was 
effected by distillation from a 12-liter flask heated internally with a 
steam coil. The vapors were condensed by passing through a !/2-in. 
copper-coil water-cooled condenser. The yield of crude pinacol 
hydrate was 52 to 58 percent of the calculated amount based on 
magnesium. About 50 kg was prepared. 

b. Pinacolone.—Pinacolone was prepared in 70-percent yield from 
pinacol hydrate by the method described in Organic Syntheses [4]. 
Double quantities were used, and the reaction was carried out in a 
5-liter three-necked flask. 

The crude pinacolone was distilled through a 100- by 2-cm column 
packed with Pyrex helices. Distillate representing 80 to 83 percent 
of the charge was collected between 105° and 106°C 

C. Pinacolyl aleohol.—A 12-liter round-bottomed flask was fitted 
with a 150- by 2-cm condenser cooled with water at 5° to 10°C. The 
flask was charged with 740 ml (6 moles) of purified pinacolone, 3,600 
ml of ether and 2,700 ml of saturated sodium carbonate solution. 
Sodium (690 g, 30 moles) cut into small pieces was added through the 
condenser as fast as possible without causing an excessively violent 
reaction. This required from 6 to 7 hours. When sodium carbonate 
crystallized out to such an extent that the mixing caused by evolution 
of hydrogen was impaired, water was added, 25 ml at a time, until 
the salt was just dissolved. 

When the reaction was complete, the layers were separated, the 
aqueous layer extracted once with 1,000 ml of ether, and the combined 
ether solutions were washed with water, dried with potassium carbon- 
ate, and distilled through a 100- by 2-cm packed column. The product 
collected boiled at 120° to 122°C. In this way 10,850 ml of pinacolyl 
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alcohol was prepared from 12,295 ml of pinacolone (88.4 percent of the 
calculated yield). 

D. Dehydration of pinacolyl alcohol—Pinacolyl alcohol was de- 
hydrated by passing it through a 100- by 2-cm horizontal Pyrex-glass 
tube electrically heated, well lagged, and filled with alumina (Hydralo, 
from J. T. Baker Chemical Co.). The alcohol was added at the rate 
of 60 ml per hour while the tube was maintained at 297° to 305° C. 
The products, olefins and water, were separated, and the olefins dried 
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Figure 1.—Fractional distillation of dimethylbutenes. 


with anhydrous sodium sulfate. In this way 7,742 ml of mixed 
olefins was obtained from 7,795 ml of the alcohol. The product was 
water-white and contained no alcohol. The mixture of olefins was 
distilled through the 150-plate bubble-cap column. One charge, the 
distillation curve of which is shown in figure 1, was separated into 57 
fractions, 27 of which boiled between 41.23° and 41.24° C (head 
temperatures), 21 between 55.56° and 55.62° C, and 5 between 73.17° 
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and 73.18° C. Not including the holdup (210 ml, 3% fractions) and 
four intermediate fractions, this represents 51.0 percent of 3,3-dimeth- 
ylbutene-1, 39.6 percent of 2,3-dimethylbutene-1, and 9.4 percent 
of 2,3-dimethylbutene-2. 

Since approximately equal portions of the three olefins were desired, 
excess 3,3-dimethylbutene-1 and intermediate fractions were passed 
through the catalyst tube previously described, filled with aluminum 
sulfate maintained’at.280° to'290° C. The first’portion of the effluent 
contained hydrogen sulfide and sulfur, but the remainder was clean. 
All olefins from the reaction were washed with sodium plumbite and 
shaken with mercury. From 3,570 ml charged, there was obtained 
3.390 ml of mixed olefins. The relative proportions, as determined 
by distillation analysis, were: 3,3-dimethylbutene-1, 7 percent; 2,3- 
dimethylbutene-1, 28 percent; and 2,3-dimethylbutene-2, 65 percent. 
These results substantiate those found by Cramer and Glasebrook [5]. 

Appropriate fractions from several runs were combined and redis- 
tilled, and samples of the middle portion of each distillate were used 
in the determination of physical constants. 


2. PREPARATION OF 2,2-DIMETHYLBUTANE 
Flow diagram: 


H, 
3,3-dimethylbutene-1 —--————> 2,2-dimethylbutane 
Raney Ni 
catalyst 


Fifteen hundred milliliters of 3,3-dimethylbutene-1 was hydro- 
genated at 1,500 lb/in.” pressure, at 150° C, in the presence of Raney 
nickel catalyst. After washing, drying, and refluxing over sodium, 
1,485 ml was carefully fractionated in the 600- by 4-cm packed column. 


3. PREPARATION OF 2,3-DIMETHYLBUTANE 


Flow diagram: 


H; 
2,3-dimethylbutene-2 —-—————> 2,3-dimethylbutane 
Raney Ni 
catalyst 


One thousand milliliters of constant-boiling 2,3-dimethylbutene-2 
was hydrogenated at 1,800 lb/in.? pressure, at 130° C, in the presence 
of Raney nickel catalyst. The product was washed, dried, refluxed 
over sodium, and distilled in the manner already described. The 
product, 987 ml, was carefully fractionated in the 600- by 4-cm 
packed column. 


4, PREPARATION OF 2,3,3-TRIMETHYLBUTENE-1 AND 2,2,3-TRI- 
METHYLBUTANE 
Flow diagram: 
acetone 
tert-butyl chloride 


ether 


> 


tert-butyl magnesium chloride 





Mg 
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P dehydrate ‘ 
2,3,3-trimethylbutanol-2 2,3,3-trimethylbutene-1 


H, 
—_————> 2,2,3-trimethylbutane 
Raney Ni 
catalyst 


2,3,3-Trimethylbutanol-2 was prepared by the action of acetone 
on tert-butyl magnesium chloride in dry ether in yields up to 30 per- 
cent. The product was recovered, dehydrated to 2,3,3-trimethyl- 
butene-1, and hydrogenated to 2,2,3-trimethylbutane in a manner 
essentially the same as that described hereinafter for the preparation 
of 3-ethylpentane. The carbinol boiled at 129° to 133° C, the olefin 
at 76° to 78° C, and the paraffin at 80° to 82°C in the preliminary 
separations. ‘ 

A portion of the 2,3,3-trimethylbutene-1 was purified by distilla- 
tion through the 150-plate bubble-cap column. From 3,000 ml of 
crude olefin 2,600 ml o 2.3,3-trimethylbutene-1 was recovered, which 
was constant-boiling (within 0.01° C) and which had a constant 
refractive index. 

All 2,2,3-trimethylbutane was carefully washed with sulfuric acid, 
water, sodium bicarbonate solution, dried over calcium chloride, and 
distilled from sodium. The distillate boiling at 80° to 82° C was frac- 
tionated in the 600- by 4-cm packed column. Six liters of material was 
prepared in this way, fractions of which differed in boiling point by 
no more than 0.013° C and which had refractive indices (n%) of 
1.3892 to 1.3894. 


5. PREPARATION OF 3-ETHYLPENTANE 
Flow diagram: 
ethylpropionate 


Ethyl bromide 


ether 





> ethyl magnesium bromide 
Mg 








2 
complex > triethylearbinol » 3-ethylpentene-2 
HC 


He 
> 3-ethylpentane 





Raney Ni 
catalyst 


Triethylcarbinol was made in 80- to 85-percent yield by the action 
of ethyl propionate on ethyl magnesium bromide in dry ether. The 
carbinol was isolated from the reaction mixture by decomposing the 
Grignard complex with ice and dilute hydrochloric acid, and extract- 
ing the aqueous layer with ether. The ether solution was distilled 
through ‘s 100- by 2-cm packed column. The carbinol boiled at 140° 
to 142°C. 
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The neseapatene-S was formed in a 94-percent yield when 
triethylearbinol was slowly distilled in the presence of 0.5 percent of 
its weight of iodine. The product was dried with calcium chloride, 
refluxed over sodium, and distilled in the 100- by 2-cem column. The 
fraction which boiled at 95° to 97° C was collected for hydrogenation. 

The 3-ethylpentane was formed in a 97-percent yield when 3-ethyl- 
pentene-2 was hydrogenated at 120° C and 1,800 to 2,200 lb/in.? pres- 
sure in the presence of Raney nickel catalyst. The product was 
thoroughly washed with concentrated sulfuric acid, water, 5-percent 
aqueous sodium bicarbonate, and finally washed several times with 
water. It was dried with calcium chloride, refluxed over sodium, and 
distilled through the 100- by 2-cm. column. The fraction which 
boiled between 92° and 94.5° C was collected for further purification. 
A total of 2,650 ml was prepared in this way. 

The material was further purified by careful fractionation in the 
600- by 4-em packed column. LEach fraction was withdrawn after 
the temperature at the head of the column (corrected to 760 mm Hg) 
had become constant. In this way 1,750 ml of hydrocarbon was 
collected at still head temperatures within a range of 0.012°C. All 
fractions of this portion had a refractive index (n7°) of 1.3928. 


6. ISOLATION OF 2,3,4-TRIMETHYLPENTANE 


This octane was prepared by fractionation of a mixture of isomeric 
octanes supplied through the kindness of Gustav Egloff, of the Uni- 
versal Oil Products Co. The isomeric octanes were first distilled in a 
batch still having a 240- by 4-cm column packed with locket chain. 
Fractions of about 4.5 liters so obtained were distilled in a 150-plate 
— still, recombined according to refractive index, and redis- 
stilled. 


7, PURIFICATION OF n-HEPTANE AND 2,2,4-TRIMETHYLPENTANE 


Samples from the purest batches of certified n-heptane (California 
Chemical Co.) and “isooctane’’ (2,2,4-trimethylpentane) (Réhm & 
Haas Co., Inc.) were carefully fractionated in a 150-plate bubble-cap 
column. 


IV. DETERMINATION OF PHYSICAL PROPERTIES 


The apparatus and technique employed in the measurements of 
physical properties have already been described in the first report 
1] of this series. 

In previous work [1, 6, 7] it has been found that data obtained in 
the distillation of a sample through the 10-in. jacketed column used 
in determining the boiling point of the sample furnished a rough index 
of its purity. The distillation data on each sample were obtained at 
two or more total pressures differing by 7 to 28 mm Hg, one of these 
values being close to 760 mm Hg. From the rate of change of tem- 
perature with pressure so found, the data obtained when the pressure 
was near 760 mm Hg were corrected to that value. As pressure 
readings were taken only at intervals during a distillation, minor 
fluctuations not recorded caused an unevenness in the corrected data 
amounting to three or four thousandths of a degree. The corrected 
distillation data obtained in the present work are plotted in figures 
2 to 9 and 11 to 13. 
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The two lower-boiling hexenes, 3,3-dimethylbutene-1 and 23. 
dimethylbutene-1, figures 2 and 3, appear to be of fairly high purity, 
By comparison, the 2 o- dimethyl butene-2, figure 4, is of very low 


purity. 
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FiaurEe 2.—Distillation of 3,3-dimethylbutene-1. 
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Fiaure 3.—Distillation of 2,3-dimethylbutene-1. 
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Fraure 4.—Distillation of 2,3-dimethylbutene-2. 
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Figure 5.—Distillation of 2,2-dimethylbutane. 


The two dimethylbutanes, figures 5 and 6, are of about as high 
purity as the lower-boiling hexenes. The 3-ethylpentane, figure 7, 
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is of very high purity, no rise of temperature being observed during 
distillation from 6 percent to 90 percent recovered. 

The 2,3,3-trimethylbutene-1, figure 8, appears to be of lower 
purity than the hexenes. Its paraffin homologue, 2,2,3-trimethyl- 
butane, figure 9, like the 3-ethylpentane, shows no rise of temperature 
in distillation above 3 percent recovered. A freezing-point determina- 
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Fiaure 6.—Distillation of 2,3-dimethylbutane. 
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Figure 7.—Distillation of 3-ethylpentane. 
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Figure 8.—Distillation of 2,3,3-trimethylbutene-1. 
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Figure 9.—Distillation of 2,2,3-trimethylbutane. 


tion on this material is plotted in figure 10. A freezing-point deter- 
mination on another sample, to which 2 percent of normal heptane 
had been added, showed the lowering to be 2.17° C per mole-percent 
of impurity. ‘ 
The distillation curves of normal heptane, figure 11, and of 2,2,4-tri- 
methylpentane, figure 12, offer no especial points of interest. Figure 
13, the distillation curve of 2,3,4-trimethylpentane, is of interest in 
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FicurE 10.—Freezing curve of 2,2,3-trimethylbutane. 
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Figure 11.—Distillation of n-heptane. 
































=. >—< 











| | | tt | 



































30 40 so 60 70 80 
RECOVERY, PERCENT 


Figure 12.—Distillation of 2,2,4-trimethylpentane. 
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Figure 13.—Distillation of 2,3,4-trimethylpentane. 
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that this material was isolated by fractional distillation from a syn- 
thetic crude containing 5 to 10 percent of hexanes and heptanes, the 
balance being octanes including the four trimethylpentanes, at least 
two dimethylhexanes, and small amounts of unidentified octanes. 
Without recourse to fractional crystallization, a product was ob- 
tained, 80 percent of which distilled within a range of 0.03° C. The 
freezing and melting curves of 2,3,4-trimethylpentane are shown in 
figure 14. 
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Figure 14.—F reezing and melting curves of 2,3,4-trimethylpentane. 


TABLE 1.—Boiling-point data on 2,2,3-trimethylbutane 





— | Observed | Computed Difference, | 
| ————_ | boiling | — boiling observed — 
P point point computed 


mm Hg “¢ °C 

740. 74 } 80.006 80. 005 
750. 20 | 80. 434 80. 436 
758. 67 | *80.819 | 80. 819 
764. 41 | 81.082 81.076 
767. 44 81. 207 81. 211 


| | | 
* This value was obtained on a sample from a different batch of 2,2,3-trimethylbutane. 














TABLE 2.—Empirical constants in vapor pressure equations 





Compound j A 
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3,3-Dimethylbutene-L. | 17.04606 
2,3-Dimethylbutene-L PaSctaccwosauseces=| Sanu 
2,3-Dimethylbutene-2_ . ee eee 17. 63641 
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3-Eth ee ee 
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»4-Trimethylpentane. | 17.25410 
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Boiling-point determinations were made on 2,2,3-trimethylbutane 
at several barometric pressures, the values being given in the first 
two columns of table 1. Since the heat of vaporization could be 
assumed to be constant over the small range of temperatures shown 
in the second column of the table, the equation 


log p=A—B/T, (1) 


where p=pressure; 7=absolute temperature; and A, B=empirical 
constants, should fit the data satisfactorily. The constants A and PB 
were evaluated from the data by the method of least squares, and 
are given in table 2 for nine of these compounds and for 2,2,3-trimethy]- 
pentane, the properties of which have been described earlier [1]. 
Values are not given for n-heptane and 2,2,4-trimethylpentane, as the 
cubic equations of Smith and Matheson [9] are more accurate. 

The third column of table 1 gives the boiling points computed from 
the formula so found, and the last column shows the difference be- 
tween observed and computed values. These differences are of the 
magnitude to be expected if the standard deviation of a boiling-point 
determination is +0.003° C. This value is only one-third as large as 
that found in earlier work with this apparatus [7]. Refinements in 
technique are believed to be responsible for the reduction. 


TABLE 3.—Properties of preparations * 





| | 

Boil- | Refractive | : 
Freez- | ing | index Density 
Compound | Tia | ge] aur | 
point | 





mm | 9% 
Hg “D 





deg/mm 
wd 4 7 Hg per deg ml deg 
~~ —115. 53) 41. 3 . 37604) 1. 37295 ’ R 30. —0. 001014 
3,3-Dimethylbutene-1_. - (rie? | tLe as, 87660) 000598| . 65198). »—. 000908 
—140.1 | 55.615). | 1. 39022)1. 
2,3-Dimethylbutene-1___!{—120 to | 
—193 55.8 | 8915 - 
2,3- Dimethylbutene-2__- f oo es eae" tne 7 “7064 \_ ae 


2,3,3-Trimethylbutene-1_ : 7 — iz | Y y ; | *Fom'| : 


| 
, , 06, 49.7311 . 368641. —..000543} . 
ee 1 eT ‘36887 =! 000541) : 


2,3-Dimethylbutane.._ _. i y = o| : “3409 

















2,2,3-Trimethylbutane._ { = -* - . 
ja 65) 93. 473 F . 39337) 1. 39084 
—118.8 | 93.3 \°], $937 

ie ty 98. 424) - 38764 /1. 38511 


—90.5 | 98.4 88774 *—. 00014} . 








{—107. 37 99.233) «(. _ 39155/1. 38915, —.000480| . 
—107. 45| 99.3 | 1.89187|......- ; 


|f —109. 97113. 391) 1. 40431|1. 40193) —. 000476) . 
L (® |us.4 |. 1 4048 | E 


2,2,4-Trimethylpentane- 
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® The values given in the upper rowbt figures for each compound are those obtained in this work and are 
believed to be accurate within a few units in the last decimal place given. The lower row of figures (in italics) 
are taken from Egloff [8] for comparison. 

» Over s wider range of temperatures. 

e At 18°C, 

4 Estimated. 

¢ Probably a typographical error. 

! “Glassy at low temperatures.” 

& Copied from reference [7]. 
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The physical properties of the compounds included in this investi- 
gation are given in table 3. The upper row of figures for each com- 
pound are the values obtained in this work. For comparison, values 
taken from Egloff [8] are give below in italics. The boiling point of 
each compound at 760 mm Hg barometric pressure was found by 
substituting in eq 1 the observed data for the compound, evaluating 
the constants A and B, and solving the resulting formula for p=760. 
The values so obtained deviate from those shown in figures 2 to 9 
and 11 to 13 by one or two thousandths of a degree, the latter having 
been corrected to 760 mm Hg by linear interpolation only. Refrac- 
tive indices and densities were obtained at temperatures close to 20° 
and 25°C and were corrected to 20° and 25° C by assuming the vari- 
ations with temperature to be linear. All weighings made in the 
density determinations were corrected to vacuum. 
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MATHEMATICAL THEORY OF IRROTATIONAL 
TRANSLATION WAVES 
By Garbis H. Keulegan and George W. Patterson 


ABSTRACT 


‘This paper is the first of a series dealing with the motion of flood waves and 
other waves of translation in open channels. The case treated is that of waves 
for which the forces of fluid friction are negligible with respect to the inertia 
and gravitational forces. The irrotational motion of a perfect liquid in a hori- 
zontal rectangular canal when the original surface is disturbed is investigated 
on the assumption that the horizontal velocity in a cross section is approximately 
uniform. The results are also applicable to motion in a canal of uniform slope 
containing water originally moving with a uniform velocity. Special emphasis 
is laid on disturbances which are propagated without change of form, and in 
these cases formulas are derived for the wave profile and velocity of propagation. 
Formulas are also derived which give the deformation, energy, motion of the 
center of gravity, and moment of instability of an arbitrary intumescence. Con- 
sideration is given to the maximum height of a wave of permanent form. For- 
mulas have been compared with the available experimental data. Of special 
interest is the comparison of the shape of the undulations composing the head of 
an initial sur6e with the characteristics of the cnoidal wave. 
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LIST OF SYMBOLS 


c velocity of propagation of long wave of negligible height and 
curvature. 
c¢, velocity of wave approaching a discontinuity in a channel] 
cross section. Also a constant of integration. 
c, velocity of transmitted wave after passing a discontinuity in 
a channel cross section. Also a constant of integration. 
cn (2, k) Jacobian elliptic function, cosine amplitude. 
total energy of a wave per unit width of channel [MLT~?}. 
complete elliptic integral of the second kind. 
with or without subscripts—functional symbols. 
incomplete elliptic integral of the first kind. 
complete elliptic integral of the first kind. 
acceleration of gravity. 
functional symbol. 
height of main body of discharge wave (section V—17) ; vertica! 
displacement of water surface. 
vertical displacement of surface as wave approaches discon- 
tinuity in channel cross section (section IV—5). 
vertical displacement of surface after wave passes discon- 
tinuity in channel cross section (section IV—5). 
maximum height of solitary wave, maximum depth of nega- 
tive wave, maximum height of initial discharge wave. 
maximum height and depth, respectively, of cnoidal wave. 
= Hy. 
maximum height of initial discharge wave. 
h’ average height of discharge wave. 
H_ undisturbed depth of liquid in channel. 
H( )_ functional symbol. 
k modulus of elliptic integral. 
l length of section of channel (section IV—4). 
l, m,n direction cosines. 
M, moment of instability [L]. 
M’;, Me, M’2, Mz racemes er unit width [MT~']. 
p pressure [ML~!T~-?]. 
Pa atmospheric pressure. 
q_ speed of particle (section II-5, 7, 8). 
q discharge per unit width in a cross section [L?T~'] (section 
V-2). 
Q total volume per unit width of a solitary wave [L*]. 
s length of are. 
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sn (z, k) Jacobian elliptic function, sine amplitude. 
t time. 
u,v, w 2x component, y component, z component, respectively, of 
velocity. 
uo 2x component of velocity at bottom of channel, also first 
approximation to velocity in cross section. 
Uo, Yo, Wo Velocity components prior to impulse (section II-6). 
u,, Ue velocities of particles in a wave before and after passing a 
discontinuity in a channel cross section (section IV—5). 
uy, U2, U1, Uz local and mean velocities in two cross sections (section V—16). 
u, velocity of particle at apex of wave (section V—19). 
U_ mean velocity of particles in cross section. 
Up mean velocity in a section where h=0. 
V_ volume [L’] (section II-3); volume per unit width [L?*] (sec- 
tions [V-3, V-—5). 
zx Cartesian coordinate, specifically, parallel to channel in 
bottom plane. 
y Cartesian coordinate, specifically, lateral to channel in the 
horizontal plane. 
z Cartesian coordinate, specifically, drawn upward with origin 
at bottom of channel. 
a=3/H?® (section V—-10); a ratio (section V—14). 
V1, Y2, ¥3 rates of angular dilatation [T-1}. 
T circulation (L?T-"}. 
€1, €2, €; rates of linear dilatation [T—}]. 
e asmall quantity. 
=,7, ¢ components of vorticity [T-!] (sections II-1, 4, 5). 
¢ dimensionless x coordinate of free surface. 
n,& 2 and zg coordinates, respectively, of center of gravity of 
volume above undisturbed surface (sections V—5, 6, 12, 17). 
dimensionless z coordinate of free surface=h/H (sections V-8, 
9); Boussinesq number (section V—16). 
three roots of the differential equation of the permanent wave. 
maximum relative height of a solitary wave (section V-—8); 
maximum relative height of a cnoidal wave (section V-9, 
18). 


maximum relative depression of a cnoidal wave. 

functional symbol. 

wave length of a cnoidal wave (sections V-9, 18). 

density of liquid [ML-%]. 

partial volume per unit width of an intumescence [L?]. 

volume per unit width of an intumescence. 

velocity potential [L?T-']. 

velocity potential prior to application of impulse (section 
II-6); velocity potential, first approximation. 

functional symbol. 

velocity of propagation of an element of volume of an in- 
tumescence. 

velocity of propagation of an element of wave height, negli- 
gible surface curvature. 

velocity of propagation of an element of wave height, appre- 
ciable surface curvature. 

velocity of propagation of a cnoidal wave. 

velocity of propagation of discharge. 

impulsive pressure [ML~-!T-']. 

gravitational potential [L?T-?]. 


I. INTRODUCTION 


The growing importance of ae pres the occurrence of floods and 


the rate of travel and height of flood waves as they pass down rivers 
has led in recent years to a marked increase in the literature on this 
subject. However, in no language can there be found a comprehen- 
sive presentation of the mathematical treatment of the problem of 
such waves. Instead, various authors at different times and in dif- 
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ferent countries have treated separate aspects of the problem, so that 
the student of the subject must read many papers in several languages 
to obtain a satisfactory grasp of the problem. There is definite need 
for a systematic treatment of the results of these different methods 
of attack, presenting in coordinated form the results that have been 
attained up to the present time. 

It is for this reason that the National Bureau of Standards, at the 
suggestion of the United States Weather Bureau, has commenced the 
preparation of a series of papers dealing with the mathematical theory 
of flood waves and other waves of translation. It is not the purpose 
of these papers to furnish a practical method of predicting the rate 
of flood-wave travel and the rate of attenuation of the wave in an 
actual river channel with its complex flow conditions. Instead, the 
purpose is to furnish a sound mathematical theory on which attempts 
to solve the practical problems of flood prediction can be based. 

The present paper, the first of the series, deals with the motion of 
translation waves in channels of uniform, rectangular cross section 
when the forces of fluid friction are negligible with respect to the 
inertia and gravitational forces. Later papers in the series as now 
planned are in various stages of completion and will deal with the 
following topics: The effect of turbulence and channel slope and con- 
figuration on the motion of translation waves; the theory of quasi- 
permanent regime and the methods of prediction of flood waves; and 
recent advances in the problem of the deformation of an intumescence. 

The different regimes of the flow of water in canals and natural 
watercourses may be characterized as uniform, varied, or unsteady. 

In the uniform regime the depth and mean velocity are invariable 
with respect to both time and space. This regime is attained in 
sufficiently long stretches of natural courses, provided the slope and 
discharge are constant, and the channel has a uniform cross section. 

In varied flow the depth and mean velocity are unchanging with 
respect to time, but not with respect to space. Examples are flow 
in the backwater region of a dam, flow in a canal connecting reser- 
voirs at different elevations, and flow in transition sections. 

On the other hand, the unsteady regime is characterized by a depth 
and velocity which vary with respect to both time and space. 
Gravity waves and capillary waves belong to this class of phenomena, 
the former being the more important of the two, considered from the 
point of view of applied science. Gravity waves fall into two classes— 
surface waves and translation waves. In surface waves, which occur in 
relatively deep waters, the motion decreases rapidly from the surface 
downward. The wave motion is sinusoidal, the individual particles 
periodically traverse closed paths, and the vertical accelerations are 
always appreciable. In translation waves, on the other hand, the 
motion of all the particles in a cross section is sensibly the same at 
any given instant. As a translation wave passes, the liquid particles 
are permanently displaced by an amount proportional to the volume 
of the wave. Waves of this type are frequently encountered in 
channels and rivers. Tidal waves are examples of the same phenom- 
enon in seas and oceans. 

Perhaps the simplest classification of translation waves is that 
obtained by considering the physical condition governing the motion. 
Two extreme cases are immediately apparent. First, we have the 
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case in Which the effect of friction is negligible. In these waves the 
volume and the length are small; and the latter fact explains why the 
effect of friction can be neglected, so that the weight and inertia 
become the factors controlling the motion. If the curvature of the 
wave surface is negligible, the velocity of propagation of the wave 
height, wo, is given by the Lagrange velocity law, w=+gH, where g 
is the acceleration of gravity and H is the channel depth. On the 
other hand, if the maximum curvature and the height above the 
undisturbed level of the free surface are appreciable, then the velocity 
of propagation of the wave height w, equals wo times a factor slightly 
different from unity, the exact value of the factor being a function 
of the height and the derivatives of the wave surface. We call such 
disturbances potential intumescences or irrotational translation waves. 
They can be generated by suddenly immersing a solid body in the 
liquid or suddenly adding a small mass of liquid; they may also be 
formed by the sudden motion of a body immersed in the liquid. In 
these cases the disturbed position of the surface is entirely above the 
original surface; that is, the waves have a positive elevation. The 
solitary wave is a familiar example. Waves can also be produced by 
abruptly removing from the channel! a small quantity of water or an 
object already in the channel. These waves have a negative ele- 
vation. Surges produced in channels by an abrupt change of flow at 
either end also belong in this category. 

At the other extreme, we find the quasi-permanent regime of rivers, 
or what are commonly called flood waves. These are also translation 
waves, but they are waves in which friction is the all-important factor, 
the effects of weight and inertia being negligible. The waves are very 
long, which accounts for the predominant effect of friction. It 1s 
illuminating to compare flood waves and potential intumescences on 
the basis of their velocities of propagation. When discussing flood 
waves, it is convenient to utilize the idea of the velocity of propagation 
of discharge, w,. According to the Kleitz-Seddon law [1]! this velocity 
is given by dq/dH, the slope of the rating curve, g, being the discharge 
per unit width occurring at the normal stage, H. In a rectangular 
channel w,=3U/2, if the Chezy resistance law is used. A comparison 
of the two wave velocities gives wo?/w,?=(4/9)(gH/U*). Since in 
rivers the ratio gH/U? is a large number, the conclusion is that wp is 
many times larger than w,. This difference of velocities is significant. 

It should be understood that in reality there exists a continuous 
transition between the two extreme types of waves described above. 
There are other phenomena occupying an intermediate position, in 
which the effects of inertia, height, and friction are of the same order 
of magnitude. These phenomena are also of engineering interest, par- 
ticularly in connection with hydroelectric plants. When the turbines 
of such plants are suddenly shut down, a positive intumescence is 
generated in the supply canal and travels upstream, increasing in 
length as it goes. A negative intumescence is simultaneously pro- 
duced in the discharge canal, and this travels downstream, decreasing 
inlength. These waves may be of considerable height (1 meter or so) 
and are attenuated by the forces of friction. Friction cannot be 
neglected, because of the long distances which the waves travel. 


epee 
' Figures in brackets indicate the literature references at the end of this paper. 
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Similarly, when the turbines are suddenly started, the same phenome. 
non occurs, but the signs of the intumescences are reversed. If the 
pressure and velocity variations caused by these waves are large, 
these effects must be considered in designing channel walls, bottom, 
and embankments. Waves like the above, except that they are 
generally propagated in still water, occur in canal reaches adjacent 
to locks when the latter are operated. The positive and negative 
waves that are generated are quite long, traverse large distances, and, 
although they undergo a deformation of form owing to the action of 
friction, the wave length remains the same [2]. There are also a few 
special cases of flood waves occurring in mountainous regions, in 
which the inertia effects are important. Massé [3] has given the 
name of average intumescences to the phenomena described above. 

The present paper deals exclusively with potential intumescences, 
Since friction is negligible for these waves, and since they are generated 
in a liquid either initially at rest or in uniform motion, then by the 
Lagrange-Cauchy theorem the motion of the liquid particles is irro- 
7 aos This fact justifies the mathematical approach selected for 
this paper. 

A ref presentation of the elements of irrotational theory is given 
in section II, primarily to emphasize the basic dynamical and kine- 
matical principles involved. Section III is a formulation of the basic 
differential equations whose solutions constitute the solution of the 
problem. 

The general theory is applied in sections IV and V, which deal with 
solutions of the equations to a first and to a second approximation. 
In the first approximation, the height and curvature of the waves are 
assumed negligible, and a complete solution is found in which the 
waves are propagated without change of form. In the treatment of 
the second approximation, which constitutes the main body of the 
paper, more general assumptions are made, leading to involved differ- 
ential equations which can be solved easily only for particular cases. 
The deformation of the wave profile during propagation and the motion 
of the center of gravity of the intumescence play an important role 
in this analysis. These points are discussed after the introduction 
of the concept of velocity of propagation of a volume element, origi- 
nally due to Boussinesq. There are three integrals invariant with 
respect to time associated with intumescences of appreciable curva- 
ture—the volume, the energy, and the moment of instability—which 
are discussed in considerable detail. 

The main emphasis is on the intumescence whose wave profile is 
propagated without change of form, and the particular solution is 
obtained for this case. If the liquid is undisturbed at infinity, we 
have the well-known solitary wave which has been investigated 
experimentally. If the liquid is in motion at infinity and this motion 
is indeterminate, we have the less familiar cnoidal waves. In the 
former case, there is an interesting principle discovered by Bous- 
sinesq—namely, for all intumescences of equal energy, the solitary 
wave is unique in possessing the minimum moment of instability. 
This theorem is significant in that it explains the frequent occurrence 
and ease of formation of the solitary wave. 

In discussing the velocity of propagation of surges, an attempt 1s 
made to determine the effect of the channel velocity distribution. 
The characteristics of the initial wave in a surge are not without 
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interest, and this wave is considered in the light of both theory and 
experiment. _ 

The concluding section touches on the problem of limiting wave 
heights and the breaking of waves. 

In the preparation of this paper, extensive use has been made of 
Lamb’s comprehensive text [4] and of a classic paper by Boussinesgq [5]. 
Whenever appropriate experimental data have been available, the 
authors have attempted to discuss the agreement between theory and 
experiment. The principal sources of data have been the well-known 
a of Bazin [6] and Russell [7, 8] and a recent paper by 
Favre [2]. 

II. IRROTATIONAL MOTION 


1. DYNAMICAL EQUATIONS OF MOTION 


We consider an elementary parallelepiped ézsyéz having its corner 
nearest the origin of the rectangular axes of reference at the point 
(z, y, 2). Designating the component velocities at the point (z, y, 2) 
by u, v, w, and the density of the liquid by p, we have for the kinetic 
reaction in the z direction of the liquid in the parallelepiped 


d 
py oubysz, (1) 


which is equal to the z component of all the forces acting on the liquid 
in the parallelepiped. If the pressure at the point (z, y, z) is p, we 
have for the z component of the force due to pressure, 


_ SP sciybz, (2) 


Denoting the potential of gravity at (x2, y, 2) by 2, we have for the z 
component of force due to the weight of liquid, 


02 
= ps ,dubyoz. (3) 


Equating eq 1 to eq 2 plus eq 3, dividing by péxdyéz, and applying the 
same reasoning in the remaining two directions, we obtain 


and 





To introduce the local variations of the velocity components, we 
consider uw, v, w as functions of z, y, z, and ¢t. With this understanding 
we write the identity 


du__Oudz, oudy , Owdz, ou 
di dz dit dy dit Os dit oF (5) 


where dz/dt, for example, represents the time variation of the position z 
of the liquid parallelepiped; that is, the velocity component wu, and 
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Ou/dt represents the time variation of wu at the fixed point (z, y, 2), 
Interpreting the quantities dy/dt and dz/dt in the same manner, we 


write from eq 5: P 
u_du, du, du, du 
dt Ot "dz? "dy * dz 
Similarly, 
dv Ow, WwW, Ov Ov 
dt dt* “ort dy ta2” 


and 





dw Ow, Ow, Ow, ow 
dt at + “dz t“dyt az | 
It should be noted that du/dt represents the acceleration of a particular 


element of the fluid followed in its motion relative to the axes 2, y, z. 
The symbol d/dt has in general this significance. 


2. DISPLACEMENT AND DEFORMATION OF A LIQUID ELEMENT 


We consider the same elementary parallelepiped, 6x 6y 5z, having 
its corner nearest the origin at the point (z, y, z). Designating the 
component velocities at some instant at the point (z, y, z) by u, », w, 
we find for the relative velocities at the point (z+ éz, y+6y, 2+62) at 
the same instant, by Taylor’s theorem, 


6 _ ou, i ou. | ou, 
u=—= ar Lr oy oz 5) 


dv 


Ov Ov 
bvu= et ay 520” 


oe ow. , oO 
Ow, wW 
a et ae el 





If we write 


ou Ov ow 
a=>’ a= Oy’ a=5,’ 


ow , Ov Ou , OW Ov , Ou 
UM=Dyt dz’ det de’ det dy’ 


_2w_d_du_dw ,_d_du 
a= Oy oz’ 7 dz or’? Oe oy” | 








the relative velocities may be written 


1 1 
bu= edt to (rsby +7262) +5 (0b2z— fey)» 


iv = edy +5 (1182+ wht) +5 (60 B82): 





1 
iw= boy (rab + yby) 5 (by — nd): 
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The following physical interpretation is now evident. The motion 
of an elementary parallepiped having the point (z, y, 2) for its corner 
nearest the origin may be considered as made up of three parts. The 
first part represents a motion of translation of the whole element with 
the component velocities u, v, w. The second part represents the 
relative displacement of the particles within the parallelepiped and 
describes its deformation. The linear dilatations «, €, €3; give the rate 
of increase of the edges of the element dz, dy, 62, respectively, during 
the motion. The angular dilatations +1, v2, v3 give the rate of change 
of the angles of the edges of the parallepiped, taken two at a time. 
The third part represents a rotation of the whole element about an 
instantaneous axis having as components the angular velocities 
(1/2) & (1/2) n, (1/2) ¢. The vector (&, 7, ¢) is referred to as the vor- 
ticity of the liquid at the point (z, y, 2). 

When, throughout a finite portion of a fluid mass, the component 
vorticities £, 7, ¢ vanish, the motion of that portion is called irrotational. 


3. KINEMATICAL RELATIONS 


Since €1, €2, €; are the rates of dilation of the edges of a cubical ele- 
ment, it can be readily shown that the change, dV, of the volume of 
the liquid mass, AV = 6z dy 52, due to the displacement of the particles 


per unit time is 
dV= (e+ +3) V. (10) 
We also have the relation, 
poV=(pot+Ap)(V+dV), (11) 


derived from the principle of conservation of mass. pp is the initial 
density of liquid and Ap the change in the density during the motion 
of the element. If the liquid is incompressible, Ao=0, and we obtain 
from eq 11, dV=0, and from eq 10 


Ou , Ov 

Ont oy 
which is a kinematical relation that is satisfied in the region occupied 
by an incompressible liquid. Equation 12 is known as the equation of 
continuity for an incompressible liquid. 

At the boundary of the liquid the equation of continuity is replaced 
by a different relation. e shall consider the relation in its most 
general form. Let F (2, y, 2, t)=0 be the equation of the bounding 
surface. Let u, v, w be the velocity components of a liquid particle 
on the boundary, having the position (x, y, 2) at time f. Rt time 
t+At the position of the saetels is 


(x+udt, y+vAt, z+waAt). 
Since the particle does not leave the boundary [4, p. 7], we must have, 
because of the equation of the boundary, 
F(x+udt, y+vAt, 2-+waAt, t+ At)=0. 
When the latter expression is expanded by means of Taylor’s theorem, 


we obtain the relation 


dF _ OF, OF, OF, OF ' 
ya” let Male yt on (13) 


which is the general boundary condition. 


ow 
ut ee, (12) 





56 =~ Journal of Research of the National Bureau of Standards  [wai,s, 


Two special forms are to be noted. At a fixed boundary 
lu+mo+nw=0, (14) 


/ 


where /, m, n represent the direction cosines of the normal to the 
boundary. At a surface of discontinuity 


U(u, — U2) +-m (vr — v2) +n(w;—w2) =0, (15) 


where the suffixes are used to distinguish the velocities at the two 
sides of the surface. 


4. IRROTATIONAL MOTION 


In a large and important class of hydrodynamical phenomena, the 
vorticities £, 7, ¢ vanish—that is: 


and 





As a result of this restriction, the expressions for the velocity com- 
ponents, wu, v, w, are very much simplified. Byis why, 8 


We consider now the function 4, defined | by < the’ ‘differential 


equation i 

—do=udz-+ vdy+wdz. (17) 
According to a theorem in differential equations, the necessary and 
sufficient conditions for d@ to be an exact differential, that is, for 
d@ to be expressed as 

__O¢ og Od » 

do=- dat at 520% (18) 

are the relations in eq 16. Comparing eq 18 with eq 17, we write for 
the velocity components of’a liquid particle of zero vorticity: 


and 





The function ¢ is called a velocity potential from its analogy with 
potentials arising in some other branches of physics. As the velocity 
potential exists if and only if the flow is irrotational, the latter flow is 
also called potential flow. 

If the liquid, in addition to being in irrotational motion, is also in- 
maamuniiae, then the velocity potential ¢ satisfies Laplace’s 
equation, 


2 
$=0. (20) 


0? ¢ , °>¢ 
“Oat t Oye t 





fen Irrotational Translation Waves 

This can be seen by substituting the expressions for the velocities 
given in eq 19 in the equation of continuity, eq 12. Thus, for 
the conditions of flow assumed here, (the irrotational motion of an 
incompressible liquid) the complete determination of the velocities 
at all points reduces to the selection of the harmonic function ¢ 
satisfying the boundary conditions. 


5. INVARIANCE OF THE CIRCULATION IN A MOVING CIRCUIT 


The following theorem on irrotational motion is of great importance; 
If the motion of any portion of a perfect fluid is irrotational at any one 
moment, it will continue to be irrotational, provided the external forces 
have a potential and the density of the fluid be either constant or a function 
of pressure only. 

oe of the methods for establishing the theorem is to consider the 
circulation around a closed curve of particles and inquire into the 
variation of the circulation as the particles of the curve move into new 

ositions. By the circulation around the circuit, C, is meant the line 
integral at any time, f¢, 


= (ude-+-vdy-+wdz), (21) 


of the velocity vector taken around CQ, a closed path of particles in the 
liquid. The circulation of the circuit of the same particles at time 
t+dt is 


r+ dt, 
where 
a= 


an) i qy ude + ody + wdz). (22) 


To evaluate dI'/dt we take an element of circuit C of length 6s and con- 
sider the time-rate of variation of the quantity uéz as the particles on 
the element move into new positions. Differentiating the product 
udr, we obtain 


5 (uéx) = Sin tu Lar, (23) 
The first term in the right-hand member contains the z component 


of the acceleration of a particle on ds, given by the first of the following 
dynamical equations of motion 
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Here p is the pressure, p the density of liquid, and the potential 
of the external forces at the point under consideration. The second 
term is readily transformed into 


d dx 
Un (6x) = 05 qa 


since dz/dt is the velocity u. Thus, 


d ia 
di (uéx) = : 


and similarly 
yy oy + viv, 


and 





: 62+ wee) 


Introducing eq 26 into eq 22, we obtain 


avd dp 1 f . 
Te QE Pa 45s Q a 


where g’?=u?+v*+w*. Since Q and gq’ are single-valued, continuous 
functions of z, y, z, the last two integrals in the right-hand member 
of the equation vanish. If p is a single-valued, nonvanishing, con- 
tinuous function of p only, or a constant, then the first integral 
vanishes also. Hence, under these circumstances 


adv = 
dt 0; (28) 
that is, the circulation around a circuit moving with the liquid re- 
mains unaltered. 

By using Stokes’ theorem, the circulation, I, around circuit C, 
may be related to the vorticity of the liquid particles on a surface S 
bounded by C. Stokes’ theorem states that the line integral of any 
vector around a closed curve is equivalent to the flux of the curl of 
the vector through any surface bounded by the closed curve. Using 
this theorem, we now write 


P= fa(lét+mn+ng)ds, (29) 


where 1, m, n are the direction cosines of the normal drawn inward 
from the surface. 

It now follows from eq 28 and 29 that the property of irrotation- 
ality of a portion of the fluid is carried with the fluid. For, if initially 
§, 7, ¢ vanish in a region containing a fluid, then the circulation for any 
circuit in the same region vanishes. This is true at all subsequent 
times in the region containing the same particles, by eq 28 above, and 
the vorticities consequently remain equal to zero. 
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The restrictions under which the main theorem, expressed by eq 28, 
has been proven must be remembered. It is assumed that the 
external body forces have a potential, and that the density is either 
constant or a function of pressure only. The latter condition is 
violated, for example, in the case of convection currents generated by 
a temperature gradient. Since in the dynamical eq 24, which were 
employed in the proof, the viscosity terms are missing, another 
restriction is the condition that the effects of viscous forces are 
negligible in comparison with inertial forces. 


6. IMPULSIVE GENERATION OF MOTION 


When the boundary conditions of a mass of liquid are suddenly 
changed, a sudden alteration in the motion may take place. The 
forces introduced are of the nature of impulsive pressures, and, during 
the time interval that the impulses last, the effects of all the other 
forces are negligible. 

Let p be the density, wo, %, Wo the components of velocity imme- 
diately before, wu’, v’, w’ those immediately after the application of the 
impulse, @ the impulsive pressure at the point (x, y, z). The change 
of momentum parallel to z of an elementary rectangular parallelepiped 
having edges 6x, éy, 6z, and with its corner nearest the origin at 
(x, y, 2), 18 pdxdydz(u’—uy); the resultant impulsive pressure in the 


ie ae 
same direction 1s —Fzdubysz: Hence, 


es | 
poxdydz (u’ — Uy) = = bxdy5z, 


u' —wW=— 


Similarly 


and 





In these relations we have a physical interpretation of the velocity 
potential, @. For, consider the case where Uo, %, Wo are zero. A 
state of motion for which a velocity potential exists could be produced 
from rest by the application of a properly chosen system of impulsive 
pressures. If ¢ be the velocity potential, the desired system of impul- 
sive pressures is@=p¢+C. The arbitrary constant in this expression 
shows that a uniform pressure throughout a liquid mass has no effect 
on the motion. If the initial state of motion be uniform, that is, 
if for the initial state of motion the liquid possesses the velocity 
potential ¢o>=wupr-+uy+ wz, then the new state of motion possessing 
the velocity ookantial ¢@ may be produced from the initial motion by 
the application of the impulsive pressures @=p(¢— %)+C. 
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7. PRESSURE EQUATION 


In a case of irrotational motion the velocity potential, ¢, determines 
the velocity components wu, v, w, but the hydrodynamical solution of 
the flow is not complete as long as the pressures are unknown. The 
determination of the pressures is made by means of the pressure 
equation, which is a dynamic relation to be derived presently. 

Since the motion is irrotational, the velocities wu, v, win the dynamical 


equations, 
du__10p_ o 


dt pox Ox 
do _ 

dt 
dw__10p_ 
dt 





where the operator d/dt is 


d oO rs) a) o , 
ait et? ayt” oe? (31) 


may be expressed in terms of the velocity potential, ¢. The dynamical 
equations are transformed into 
oo ,10,, 2 1 Op 


~Sebl 3 Or! Bap Oe’ 


2 


—“— Sy poy 
10 ,___ 02 1 op. 
~ 2 p 02) 


¢,10 
—dyoit 2 dy! 


Oo ¢ 
— S20 2 d2! 





Multiplying the first by dz, the second by dy, the third by dz, adding, 
and then integrating the sum, we obtain the integral 


Ss 2— (1/2) +f), (33) 


where {o is an arbitrary function of time. This equation makes it 
ossible to determine the pressure, once the velocity potential ¢ is 
own; it is called the pressure equation. 
For an incompressible liquid, 


B_8—9—S0 +f). (34) 


To determine the arbitrary function f(t) we must know p and ¢ as 
functions of t at some definite point in the region occupied by the 
liquid. Generally this point is one where 0¢/0¢ vanishes, q is a con- 
stant, and p is the atmospheric pressure pag. 





i 


oe SS ee 


=.>0—lC SSS: 
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8. SOLUTION OF PROBLEMS IN IRROTATIONAL MOTION 


The mathematical procedure for solving a problem of the irrota- 
tional motion of an incompressible liquid bounded in part by fixed 
boundaries and in part by a free surface consists in finding a har- 
monic function, $(z,y,z2,t), which satisfies Laplace’s equation (eq 20) 
in the region occupied by the liquid, and the following surface condi- 
tions. 

At the free surface, if it exists, there are two conditions to be satis- 
fied. If F(z,y,z,t) is the equation of the surface, we have 


OF, OF, OF, OF __ 
n° Se (13) 


where wu, ¥, w are the velocities at the surface. The second condition 
states that the pressure is atmospheric at the free surface; we then 
have from eq 34 


pelo= 80-4 +f(t). (35) 


~ 


At the fixed boundaries, 


Od _ 
55=0, (36) 


n being the normal drawn to the boundary. 

The function, ¢, which satisfies eq 13, 20, 35, and 36 is unique, 
except for an additive constant. For references to proofs see {4, 
p. 59]. 


III. FORMULATION OF THE WAVE PROBLEM—WAVES IN 
WATER OF SMALL DEPTH 


We consider a rectangular horizontal channel of uniform cross- 
section and infinite length, containing a perfect liquid initially at rest. 
We assume that some external agency acting for a short time produces 
a deformation of the free surface of the liquid in the channel. If the 
action of the external agency is such that the associated forces exerted 
on a portion of the water in the channel are essentially in the direction 
parallel to the direction of the channel and are constant throughout 
the whole depth of the liquid, the particles in a normal section under 
the deformed surface will be moving with a constant velocity. Under 
these conditions the deformed free surface constitutes a translation 
wave. 

The liquid of the channel being initially at rest, we infer from the 
theorems previously explained that the resulting motion of the liquid 
under the deformed free surface is irrotational and possesses a velocity 
potential, @. Let the axis of z be drawn in the bottom plane of the 
channel parallel to the channel axis; let the axis of 2 be drawn upwards. 
Denoting the velocity components in the direction of z and z by u 
and w, respectively, we have 
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od 


t=—- >= 


Ox 


Op 


w=-—-=s? 
Oz 


(38) 


where ¢ is a function of x, z, and ¢ and satisfies Laplace’s equation 


a, 2 
ot oe a 


The pressure equation (35) reduces to 


o 1 ; 
Bt —ge—x(w tw) + FO); (40) 
where g is the intensity of gravity. 

Let the z coordinate of the free surface be denoted by h+-H, where 
AZ is the initial undisturbed depth of water in the channel, and h is the 
elevation (or depression) of the wave. 

To evaluate F (t), let us consider the point at infinity on the free 
surface; at such a point 0¢/0%t=u=w=h=0, z=H, and p=pa. Sub- 
stituting in eq 40 there results 


F(t)="*+ gH 


At the surface, since p= pa, we obtain the dynamical surface condition 
from eq 40, neglecting the effect of surface tension: 


O¢ _ gh (w+-w*)=0, =H +h (41) 

Although the last equation has been derived for the case where the 
liquid at infinity is at rest, it is equally applicable to the case where 
u?+w? and 0¢/ot vanish at the points on the surface where h=0, as, 
for example, in an intumescence consisting of an infinite number o! 
undulations of the same dimensions moving without change of form. 

We obtain the kinematic surface condition, that a fluid particle of 
the free surface always remains on the free surface, by substituting 
F(z,y,2,t)=2z—H—h=0 in eq 13, obtaining 


oh, Oh 
w= st us 2z=H-+h. (42) 


At the bottom surface of the channel the velocity component normal 
to the bottom vanishes, that is, at z=0, 


oo _ 
5o=0. (43) 


In view of this last condition, we adopt for the potential ¢ the power 
series 


= Dionc", $=0, (44) 





an 


Al 
th 


an 
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where ¢, is @ function of z and ¢ only. Substituting this expression 
for ¢ in eq 39, and equating the coefficients of 2" to zero, term by term, 
we arrive at the convergent series 
2? oO? do gt ot do 2 °° ry 
o=bo— 7 Oat tH] Ox Gl Ow Ts (45) 








Differentiating with respect to z, 
oo = Ph 2d 2° 0° go 
oz oa? 3! ox 5! oz® 
The function ¢, as stated above, depends only on z and t. If we 


denote the velocity of particles of water at the bottom of the channel 
at z by Uw, then , 





(46) 


b0= ($0) »— ik oes dz= (0) e+ {- Ulex, (47) 


where (¢o)« is the value of gat z=. Since the liquid in the channel 
at infinity is at rest, we may take the value of (¢). as zero: then, 


de f "edi (48) 


In the above formulation, h and ¢p are functions specifying a trans- 
lation wave. They can be determined by integrating the two simul- 
taneous differential eq 41 and 42 representing the dynamical and 
kinematical surface conditions. 


IV. LONG WAVES OF NEGLIGIBLE HEIGHT AND 
CURVATURE 


1. VELOCITY OF PROPAGATION 


If uw? and w’ are negligible in comparison with gh (assumptions I, 
If) and wdh/dz in comparison with Oh/dt (assumption III), then the 
surface conditions given by eq 41 and 42 reduce to 


o¢ 


gh—>; 


db , Oh 
+ 5=0. (50) 


Also, if in the expansions of ¢ and _0¢/0z given in eq 45 and 46, only 
the first term is retained (assumption IV), we have 


$= do, (51) 


and 
Ob 2? do 
2 Oat _ 
198881—40——-5 
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Here, then, we have the basis for a first approximation to a solution of 
a wave configuration. The horizontal component of the velocity iy 
a normal section is constant and equals the bottom velocity 1,. 
the vertical component varies linearly with 2z. 
Substituting eq 51 and 52 in eq 49 and 50, we obtain 


and 


0? dy) , Oh __ 
—H Or + 3=9 


assuming that / is small in comparison with H and hence can be 
neglected (assumption V). Eliminating ¢ from eq 53 and 54, we 
obtain the differential equation 

(55) 
where 

(56) 
Eliminating / similarly and writing 


____ Odp 
U=— ’ 
there results for uw the relation 
0? u role ™ 
a og (58) 


Equations 55 and 58, which are well known, possess complete 
solutions, which also satisfy eq 53 and 54: 


h 


Aa~h@e—et)+hetret), (59) 


and 


= i(x—ct) — F,(x+et), (60) 


where F; and F, are two arbitrary functions. 

These results have a very simple geometrical interpretation. 
Consider the motion represented by the function F, in eq 60. When 
z and ¢ are increased by cr and 1, respectively, the value of the func- 
tion is not changed; the disturbance A which existed at the place z 
and at time ¢ has moved during the time interval 7 to the place 
z+cr. Therefore, the function F; represents a wave traveling in the 
positive z-direction without change of form at a constant velocity 
of propagation? c=./gH. 

Similarly, F, represents another wave propagated in the opposite 
direction with the same velocity. Any motion whatever subject to 
the conditions assumed above for the approximate solution will be 
made up of waves of these types. F; and F, are determinate when 


* The question of the meaning of velocity of propagation when the wave does not retain its form will be 
discussed later in section V-1. 
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; and Uw are known for every value of ¢, and at a single point z=2,; 
o when A and uw are known for every value of x and for a single 
time t=¢). 
If h and wo represent the surface elevation and velocity of particles 
in either the progressive or regressive wave (F; or F,=0) it is seen 
from eq 59 and 60 that 


h 
Up = C7 (61) 


where the plus sign applies to a wave moving in the positive x-direction 
(h= =0) and the minus sign to a wave moving in the negative z-direc- 


tion (F, =0). 
We shall now examine the restrictions assumed in obtaining our 


approximate solution: 
w<gh, (I) 


w<gh, (ID 


won oh 
Oh gh. ty 


The two series converge rapidly; that is, 


2? 0*o 
oT) “=s <<, 


ne. V) 


It can be demonstrated that the following physical restrictions lead 
‘0 the above mathematical assumptions: 

(A) The ratio of the wave amplitude to the depth of the channel 
issmall: (h/H)<€ 1. 

(B) The slope of the wave is small: (0h/Or) <1. 

oy -_ radius of curvature of the waves is large with respect 
to H*/2 


2. ENERGY OF AN INTUMESCENCE 


A single wave is called an intumescence. The potential energy 
per unit width of an intumescence, at time t, due to the elevation or 
lepression of the fluid surface above or below the level of the previ- 
ously undisturbed liquid in a channel is 


ge” h7dz, (62) 


where 2; and 2 denote the extremities of the intumescence. The 
kinetic energy per unit width of the liquid under the free surface of 
the intumescence is 
z 
eH ( asde, (63) 


‘ 
y zi 
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Substituting for % its value in eq 61, eq 63 becomes 
eg” 72 
2 f" h? dz. (64) 


Comparing eq 64 and 62, it is seen that the potential and kinetic 
energy of an intumescence are equal and that the total energy is 


1 
oo. h? de. (65) 


The integral in eq 65 without the multiplier gp will be of specia! 
interest, and will be subsequently referred to as the energy, E, of 
an intumescence. 


3. DISPLACEMENT OF THE PARTICLES 


The liquid particles in the cross section at z, before the arrival oj 
an intumescence, are at rest. At the moment when the intumescence 
is passing, the particles in the section are moving with the velocity 1, 
to a first approximation. Thus, the volume of flow through the 
section during the passage of the intumescence is 


ta 
V= f un,iHdt. (66) 
Since 
c dt=dz, 


we have 


V= ° wd. 


Substituting for uy from eq 61, we obtain 


V= f . hdz. 
1 


The volume of flow through a cross section equals the volume of the 
intumescence. Hence, during the passage of an intumescence, the 
particles in a cross section are moved forward a distance 6z, the mag- 
nitude of which is given by the ratio of the volume of the intumescence 
to the area of the channel cross section. 


4. EFFECT OF AN ARBITRARY INITIAL DISTURBANCE 


If for t=0 the velocity of the liquid particles'in the channel and the 
form of the disturbed surface are known: 


U 
= 9), —ocrc+o, 


=y(z), —~<2z<4+o, 





—- Irrotational Translation Waves 


then the functions F; and F; occurring in eq 59 and 60 are 


F(e) =510() + ¥@), 
and , (72) 
P(t) =5{—0(2) +0.) 


To determine the form of the wave jprofile, h/H=y/(z,t,) for any 
subsequent time, ¢=t,, the curve representing F(z) is displaced in the 
positive z-direction by a distance ct,, the curve representing F(z) by 
the same distance in the opposite direction, and the ordinates of the 
two displaced curves are added algebraically. 

Consider disturbances confined initially to a finite segment of the 
channel of length 7. Resolution of the initial disturbance into F; and 
Ff, indicates that after the lapse of the time interval r=c/2, two waves 
of equal length / will be traveling in opposite directions. In general, 
resolved waves will have unequal heights at corresponding points. 
If, however, in the initial state, u=0, that is, 6(r)=0, then F,(z)= 
F(z), and the two component waves have the same shape. If, on 
the other hand, w= -+ch/H initially, then either F, or F, is zero, and 
the result is a single wave traveling in one direction. 


5. REFLECTION AND TRANSMISSION OF WAVES 


A wave is reflected when it reaches a barrier in a channel. We take 
the origin, z=0, at the barrier and note that u=O at the barrier. 
Utilizing this fact in the general eq 60, we obtain the condition that 


F,(—et)=F,(+et) (73) 


for the incident and reflected waves. 
Applying this last condition to the general eq 59 and 60, we obtain 


f= F(ct—2)+Fr(ct+2), (74) 
and 

= F,(ct—2) — F(et-+ (75) 

rae x) —F,(ct-+2), (75) 


which are the equations for reflected waves or for identical waves 
traveling in opposite directions. Here, if F.,(ct+2) is taken as the 
incident wave, /,(ct—z) will be the reflected wave. At x=2,, and 
t=t,, the elevation or depression of the wave is 


Fa Fct+2). (76) 


The corresponding elevation or depression of the reflected wave will 
reach the point z=z, at the time t=t,+-2z,/c, and the value of the 
elevation or depression is 


A= Falelt, + 20/0} —2) =F, (ct, -+2). 
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Hence, elevations and depressions are reflected without change. |; 
can be shown in a similar manner that corresponding velocities arp 
reversed with their magnitudes unchanged. 

When a wave reaches a point in a channel where there is an abrup; 
change in the cross section of the channel, either in width or in depth, 
a part of the approaching wave is transmitted and a part is reflected, 
Taking the origin at the point of the abrupt change of section, we 
may write from the general eq 59 and 60 for the side from which the 
wave approaches, moving in the positive z-direction, 


= Fe—at) +fetts), as 


u 

oF @—at) —f(et+z), (79) 

in which the first term in the right member of each equation applies 

to the approaching wave, and the second term to the reflected wave. 
For the transmitted wave we have 


h 
H,=9 est), (80 


and 
2 9(z—et). (81) 
C2 


In the above equations H, and H, are the depths of the channel 
on the two sides of x=0, and c, and c, the velocities of propagation 
for the corresponding portions of the channel. At z=0, h,=ha, since 
the surface is continuous, and therefore, 


A,[F(— et) +f (e,t)] = Hog (— eat). (82 


Also, the conservation of matter requires that B,Hju,=B,H2u2, where 
B, and B, are the widths of the channel on the two sides of s=0. This 
condition gives 


B,Mye,|F (—e,t) —f (ext) ] = BeHoe2g (— ent). (83) 


We deduce from eq 83 and 84 the ratio of the elevations in the cor- 
responding parts of the reflected and incident waves and of the trans- 
mitted and incident waves. These are 


f Bia — By, 


PF Boe, Baer és 


9g oon 2B,c,H, (85) 
Fo Be +By¢:)H,’ 





respectively. 

A wave traveling in a channel connected to a body of water in a 
basin of infinite depth and width may be considered as a special case. 
For this case it is evident by inspection that g/F=0 and f{/F=—-1. 
Hence the incident wave is completely reflected at the juncture, and 
the elevations are reflected as depressions. 





Keulegan Trrotational Translation Waves 69 


Patterson ] 


v. WAVES OF APPRECIABLE HEIGHT AND CURVATURE 
1, TWO DEFINITIONS OF VELOCITY OF PROPAGATION 


For long waves of infinitesimal height, the concept of velocity of 
propagation entails no ambiguity, because such waves are propagated 
without change of form. However, when a traveling wave is under- 
going deformation, the meaning of velocity of propagation must be 
clearly defined. We may adopt either of two distinct definitions. 
The first definition will be developed on a purely formal basis. In a 
given wave system we consider the quantities Oh/Ot and Oh/dz for an 
arbitrary z and ¢, and write the following relation from eq 13: 


oh oh = 
Fy aes: oh h=h(a, t). (86) 


In effect, this relation defines the quantity w;, which we may call the 
velocity of propagation of a wave-height element. When the traveling 
wave is undergoing a change of form at a given instant, t, w, is a 
function of h. 

A second definition, which is due to Boussinesq, considers the vol- 
ume of an intumescence. The volume of an intumescence is the 
quantity of liquid in the entire intumescence between the free sur- 
face and the undisturbed primitive level of the liquid in channel. We 
have, therefore 


_ i hee, (87) 


where 2% is the rear of the intumescence, assumed moving in the 
positive z-direction. The intumescence may be limited posteriorly; 
that is, for all values of zS 2, h is negligible, and in this case 2 has a 
finite value. On the other hand, the intumescence may be unlimited; 
in this case we consider the volume, 2, included between the head 
and an abscissa, 2), considerably removed from the head, and moving 
in such a way that 2 remains constant. 

We now consider a partial volume, o, of , between the head of the 
intumescence and an abscissa, 2x: 


on { ” nde, 25%, (88) 


We assume that a plane originally at z moves in such a way that it 
always has in front of it the same partial volumec of the intumescence. 
We refer to the velocity at which this plane moves as the velocity of 
propagation of an element of volume of an intumescence, and denote 
this velocity by w. 

The relation between w and h may be obtained from eq 88 by 
differentiating the right-hand member of the equation with respect 
to t and equating to zero. However, it may be more instructive to 
obtain the desired relation in the following manner (see fig. 1). The 
curves abc and a’b’c’ represent the free surface of the moving intu- 
mescence at the times ¢ and t+ At, respectively. A plane at z having 
in front of it, at time ¢, the partial volume, o, of the intumescence, 
moves to the position z+-wAt at time t+ At and still has the same par- 
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> +5 





x + Wat 
Figure 1.—Diagram illustrating the velocity of propagation of a volume element. 


tial volume, o, in front of it. The areas A and B are equal to each 
other. Consequently, , 
"2 


hot} Atdr. (89) 


Dividing by At, 


ho= [ win, LS%, 


and then differentiating with respect to z, we obtain 
oh , o 
ot +5; (he) =0, (90) 


which is the relation connecting w with h. 
Eliminating 0h/dt between eq 86 and 90, we find that the difference 
between the two velocities of propagation is 


Ow /dh 
ow h$2 | (91) 


If a wave is moving without change of form, wis constant, and the two 
definitions give identical values for the two velocities of propagation, 
a fact that is self-evident. 


2. MEAN VELOCITY OF PARTICLES 


_ If U and q, respectively, denote the mean velocity of the particles 
in, and the discharge through, a cross section at z, then 
q=U(H+h) (92) 


Consider the inflow and outflow in a time interval, At, for the vol- 
ume between two cross sections separated by a distance, Az; Az 
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and At having the same order of magnitude. The mean discharge 
at z in the time At is 

At oO 

a+s 3 


omitting infinitesimals in At of the second and higher orders. The 
volume of the inflow in the time Af is then 


At)? d 
gar AO s 


Similarly, the mean discharge at z+ Az in the time At is 
oq At 0g 
It 9, “7+ 7 oF 


omitting infinitesimals in Az and At of the second and higher orders. 
The volume of the outflow at z+ Az in the time At is 


gat +2 anat+ Se o =, 


and the accumulation of liquid in the volume is 


_°9 
“— AvAf, 


The increase of volume of the liquid is also 


oh 
= AzxAt, 


omitting infinitesimals of the same orders as before. 
Therefore, 


Og oh 


7 oe 


oz st 
or from eq 92 


2tu7+iy}=—2. 


Substituting from eq 90 
r) _) 
ag VA+A=5 (ho) (94) 


which when integrated becomes 
U(H+h)=hwo+C. (95) 


Since the liquid is assumed to be undisturbed at infinity, the constant 
of integration vanishes, and the equation for the mean velocity is 


U= ey (96) 
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3. EFFECT OF CURVATURE AND HEIGHT ON VELOCITY 
OF PROPAGATION OF A VOLUME ELEMENT 


When the ratio h/H and the curvature of the wave profile are appre- 
ciable, the first approximation is no longer satisfactory, and in such 
cases recourse to a second approximation becomes necessary. The 
method which will now be presented is due essentially to Boussinesg. 

In the surface conditions, eq 41 and 42, we neglect w* and replace ;; 
by the value of uw in the first approximation, eq 61, replacing c by its 
value in eq 56: 

Uu= H- (97) 


These conditions now become 


0d 197. | 
gh— RD: Pax(), (98) 


and i 

O¢ , Oh g , Oh _ 
Retaining the first two terms in the expressions for ¢ and 0¢/02, eq 45 
and 46, and neglecting / in the terms having the smaller values, we 


obtain 


H? & ; 
d=b- > SS" (100) 


) 0? HP? 0 
Be OTD Sat G ae at 
We differentiate eq 100 with respect to t, replacing 0¢/0# in the second 


term of the right-hand member by its value in the first approximation 
given by eq 53, that is, by gh. In this manner we obtain 


0 Od H? oh ; 

Ot ot 2% or ante) 
We next consider eq 101. In the right-hand member there are two 
small terms. In these terms we replace 0¢)/0z by its first approxima- 


tion from eq 57 and 61, that is, by — +}. In this manner, re- 


stricting ourselves to waves propagated in the positive z-direction, 
we obtain ae 
Od _ __ 772" bo, 9, Oh EP /g Hh, ) 
oe ee ty Aas 6 VA oe — 


Substituting from eq 102 and 103, we obtain from eq 98 and 99: 
Odo, Mh? (7708 hk] _ 
fh— st qo Ae a ie ala 


and 
0? do — of? Meh] _ 
Hye t+ v9Had Ae Oa |=? (105) 


oh 
ot 





‘nes se 


Int 


or, 


'Th 
that tl 


from ¥ 
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To eliminate ¢) between these, we differentiate eq 104 twice with 
respect to x and multiply through by H, and differentiate eq 105 with 
respect to t. In place of actually differentiating with respect to ¢, 
however, we differentiate with respect to x and multiply? eS V9H. 
Subtracting eq 105 from eq 104, there results:_ 


iy 
oh ? 3 Phy 
yi 5H or 


which is the fundamental wave equation to the second approxima- 
tion, in which the effects of the wave height and of curvature of wave 
profile are included. 

Although this equation has been derived for the case where the 
liquid in the channel is undisturbed at infinity, it remains valid for the 
case where the liquid at infinity is in motion, provided u=¥g/H h, to a 
first approximation (eq 97). 

Equation 106 becomes integrable once with respect to z when dh/dt 
isintroduced from eq 90. Performing the integration with respect to x 
and remembering that A and its derivatives vanish for z= ©, we 


_ ol h?  H? oh 
L 3 
(lus) + 9H. + gH ost <3 )= 0; 


2 
=H 4 + 9H (106) 


(107) 


which permits the determination of w, the velocity of propagation of a 
wave-volume element, in terms of h and 07h/0z’. 
We adopt as an expression for w the product 


w=7gH(1+6), (108) 


where € is a small quantity. Multiplying the two members by A and 
liferentiating with respect to t, we obtain, remembering that « is a 
small quantity and applying the same reasoning used in differ- 
entiating eq 105, 


F (hee) = VHS — gH (ch). 


Introducing 0h/dt from eq 90, this may be written, 


© (hws) =— VgH2 (hes) —gH2.(eh); 


or, in view of eq 108, 


Shes) = — gH [(1+-26)h. (110) 


' The justification for replacing the operator 0/ot by —-¥gH 2/dz can be seen from the following: Assuming 
that the wave retains its form, that is, that dF/dt=0, where F is a function of (r—ct), we have the relation 


‘tom which it is obvious that 


Fe oe ae 


While the present discussion relates to waves that may change their form, the above relation can be applied 
‘osmall terms in the wave equation, to a first approximation. 
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The substitution of eq 110 in eq 107 makes the latter expressio, 
immediately integrable with respect to x. Integrating, putting the 
constant of integration equal to zero, since A and its derivatives 
vanish at z= ©, and then dividing by h, we find 


3h , Hoh 
2e=5 AT 3h oe (111) 


Hence 


—(., 3h, Woh 
o=VoA( 143 at oh oe ?, (119) 


which is the expression for the velocity of propagation of a volume 
element of an intumescence moving in still water. 

» We next derive the expressions for w%, the horizontal velocity com. 
ponent at the bottom, and u, w, the velocity components at a point 
having the ordinate z. To determine up we consider eq 105, replacing 
h/t by —O(hw)/dz from eq 90, and then integrating with respect to 
z, the constant of integration being equated to zero. In this manner 
we obtain, since u=— 0¢)/ dz, 


_./9(,_ 2 Perk 

w= A(h 4At 3 on) (113 
From eq 45, after differentiating with respect to z, and from eq 4f, 
retaining only the first two terms of the right members, we obtain: 


2 »2 
u=h—5 a, (114) 
and 

OU , FO Us 
eee ae 


which give, after substituting for % from eq 113, and neglecting quan- 
tities of small order: 


(115 


aan: H 2N0hA 
u=V9H\ aE t\ 3— on) on |’ 


al 4 h \oh , 1 2 \oh * 
w=—2 9H (Fr—aen oat a 22H or} ai 


4. DEFORMATION OF THE WAVE PROFILE 
The expression for w, eq 112, and the relation in eq 90 permit the 
evaluation of the change in A occurring during a short interval of time 


to a sufficient approximation. Denoting gH by wo, eq 90 may be 
written as . 
re) re) 


Se bovge + 2h(o—an)] (us) 


(116 


and 


The integral of this expression is 


h=fle—wat)—12-[h(w— a), (uy 
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when tf is small. This may be proved by differentiating eq 119 and 
substituting in eq 118. In the process of differentiating the second 
term in the right member, the relation 0/0¢=«)0/dz is used, this being 
valid if ¢ is small. Eq 119 enables us to trace the wave profile curve 
as it would be for time ¢,-++#, if it is given for time 4. 

We next take an element of the intumescence, do, having the 
height 4 and base dz and follow its motion. We denote the time-rate 
of variation of A during the motion by dh/dt. From the definition of 
» we obtain 

dh Oh, oh 
di ot “or (120) 


Multiplying by / and using eq 90, we obtain 
dh %) 


niece Sis ais Pb Ron 
dt h Or 
Introducing the value of w from eq 112, we obtain 


dh —,| 3h?0h H*oho?h , hH*d% h 
hgi=— VIF aH 6 da00*t 6 az) 


which may be written, after dividing through by h, 


dhs» 1dfh? HP %h (ah? 
Gi — Raden G Lhe! (Se) |p 2 


This expression gives the rate of variation of the height of an intumes- 
cence element. 


5. MOTION OF THE CENTER OF GRAVITY OF A WAVE 


The last section treated of the motion of an element, do, of a wave. 
We may also discuss the motion of a volume, V, between two abscissas, 
4, and 2 where 2,>2,. Then 


v=[" ao, 


Let € denote the horizontal coordinate of the center of gravity of 
the volume and 7 the elevation of the center of gravity above the level 
of the undisturbed liquid in the channel. By definition, 


Za 4 
ado f do 
= “5 , and 7= —* (123) 
Differentiating with respect to ¢, we obtain 


ts adh 
do =~ do 
d fre dn Jz, dt 
a Pe. , an a a (124) 
which are expressions for the velocity of the center of gravity of the 
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Ze Z3 
volume. The two integrals | wdo and oh a0 are readily evaluated 
a1 


Zi 


by means of eq 112 and 122. Writing hdz for do, we obtain 
as =[,,, 3 {* FT? / oh\® 
{. idem VgH| V tan, hdo+ (3 |; 


adh y  Pi.eh oh a 


— 2 
Jato — Vo ty G (hot (35) il 


When these expressions are substituted in eq 123 and 124, we obtain 
for the velocity components of the center of gravity 


dé p—f. . 39 , H2(an\ mn 
F— Voll 14+ 35+Gp + | (125) 


and 
dn ___ VgH he =| rh _ oh , 4 oe 
dt 2V\4H' 6 hs (3 ‘ A126) 


Equation 125 simplifies if the slopes 0h/dz vanish at the end points 
of the intumescence considered, or if the volume V is large, yielding 


(5;) =9(+3n). (127) 


Equation 126 likewise simplifies if applied to the entire body of an 
intumescence. If the intumescence is unlimited, V is infinitely large, 
and the right-hand member vanishes. It also vanishes if the in- 
tumescence is limited, for then both h and 0h/dz vanish at the anterior 
and posterior ends of the intumescence. Therefore, the center of 
gravity of an intumescence is maintained at a constant height above 
the primitive free surface, and the velocity of propagation of this 
center of gravity is constant. [5, p. 84]. 


6. VARIATION IN ENERGY OF AN INTUMESCENCE 


It may be inferred from the invariability of » that the integral 
Sh do, when applied to an entire intumescence, is independent of 
time; that is, the integral 

[owas 
Zo 


is independent of time. When this integral is multiplied by pg/2, the 
product represents the potential energy of the intumescence, eq 62. 
On the other hand, the kinetic energy of the entire liquid is given by 


the double integral 
© H+h 
p/2 f de , (u?-+-w?)de. 


The velocity components, w and w, are given by eq 114 and 115. 
From these equations, neglecting quantities of small order, we obtain 


H+h hk?  H*0/,0h 
[OP cetwde=ol tet geet FS (128) | 
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which we multiply by pdz/2 and integrate between the limits z=z, 
and z= @. The integral hoh/dzx vanishes at infinity and also, when 
the intumescence is limited, at 2. It also vanishes at 2 when the 
intumescence is unlimited, because this integral is negligible i in com- 
parison with the first two. terms owing to the length of the intumes- 
cence. Thus the kinetic energy of the wave system is 


a {. h? det fy arr) A’ dz. 
z ea) 
The total energy of wave, /, consequently, is 


nsliioum {. Wt d+ Ph” he de. (128) 


It consists of two terms, of which the first is independent of time, 
and the second which is small in comparison with the other, may be 
shown to vary very slowly with time. 


7. THE MOMENT OF INSTABILITY OF AN INTUMESCENCE 


It was seen that every intumescence is characterized by two 
integrals, H= fhdo and V=fdo, which are independent of time. 
(See eq 65, 69, and 88). In addition to these, there also exists another 


ol 3h3 
M,= ={" I ny rp |e (129) 


integral 


which is independent of time. The invariant character of this 
integral was discovered by Boussinesq, who referred to it as the 
moment of instability [5, p. 97]. 

We consider the — 


M,—jE= = (" (3 — — 5p fee, (130) 


where £ is the integral representing the energy of the intumescence. 
Differentiating both members of the equation, we obtain, since 


dE/dt=0, 
dM, _ sy - 3h°_ 6h?| 19 oh ah 28 
dt or) Ae HPI? * “or dt 


oh OW | th 
-{- Ot 12 pnt Ppt? So 


(131) 


The term in the brackets vanishes when the intumescence is limited ; 
in the contrary case it is negligible in comparison with the integral, 
owing to the great length of the intumescence. Therefore, 


dM, _—(*dh/,, h , OW, Oh ) 
a =-{' Si 12 pt pt? = (132) 


or, in view of eq 112, 


dM. 12 "< 
‘ — ~~ hed. 133 





78 Journal of Research of the National Bureau of Standards {vas 
Now, introducing dh/dt from eq 90, and then integrating, we obtain 
dM, __6 

at ohh 

If the intumescence is limited, the right-hand member vanishes: if it is 
unlimited, this member is a quantity of the order which is neglected 


in this analysis. Therefore, to this approximation, the moment of 
instability of an intumescence is independent of time. 


(hw) 2,2 (134) 


8. SOLITARY WAVES 


The effect of height and curvature on the deformation of waves 
was discussed in the previous section on the basis of the assumption 
that deformation actually occurs. It may well be that configurations 
exist with appreciable height and curvature which are propagated 
with no deformation whatever. As a matter of fact, such waves were 
first discovered by experiment. The solitary wave of J. Scott Russell, 
consisting of a single elevation whose height may be of the same order 
of magnitude as the depth, may travel considerable distances without 
appreciable deformation [6, p. 505; 7, p. 324; 9, p. 289]. This property 
is called longevity and the law governing it may be derived from eq 106. 

Putting new variables 


(x—wt)=5H, 
h=nH, 


Gi) 


into eq 106, we get, since A on the basis of the assumption 


ot Ox? 


that the wave does not deform, 


2? 
os 


2 
6b1— 91-251) —0, (135) 
which when integrated twice gives 
2 0? n 
6bn—9n Eee ted Ta. (136) 
Since 7 and x are finite, c,=0. If we restrict ourselves to the 
2 
case where at {= ~, n=0 and SrI=0, then c.=0 and we have 
6by—9nt 251 =0. (137) 


Multiplying by a d¢ and integrating: 





Keulegan 
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Sbat—39=( 32 A (138) 


Denoting the maximum height of the wave by /,, and the corres- 
ponding value of » by m, 
h, € 
a W (139) 


At the crest of the wave, n7=m, and the wave surface is parallel to the 
bed, dn/dg=0. Substituting these values in eq 138 we obtain the 


relation between b and m, 
ow 
n=b=(4/-1). (140) 


Replacing this value of 6 in eq 138: 


dn\ 
7) laine, 
and integrating twice: 


=m sech{ 4 Sat | 


Reverting to the original variables, we obtain the equation of the wave 
profile 


h=h soch’| fF (e—at) } (142) 


The free surface is symmetrical about the normal at z=ot, and 
approaches the primitive surface as c—wt approaches + (fig. 2). 








Ficure 2.—Diagram of a solitary wave. 


The velocity of propagation, w, which may be interpreted either as 
that of the volume element or height element since the wave does 
not deform, is obtained from eq 140: 


w=g(H+h), (143) 


which is the formula adopted by Russell [7, p, 328; 8, p. 423], on 
empirical grounds and confirmed by the experiments of Toc [6, p. 
496, 515). 


198881—39——-6 
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A comparison of theory and experiment, based upon 117 observya- 
tions by Bazin [6, p. 510-514], shows the error of the formula to be 
—0.91+5.10 percent. Similar calculations based upon 77 observa- 
tions by Russell [8, p. 440-442] give +0.78+2.30 percent. The + 
term indicates the standard deviation, and this usage will be follow ed 
in subsequent comparisons. It will be noted that in each case Zero 
error is well within the limits determined by the standard deviation. 

From eqs 127 and 143 we will now show that the center of gravity 
of a wave propagated without change of form is h,/3 above the primi- 
tive surface, where h, is the maximum height of the wave. Since the 
wave retains its form, 


a, (144) 


and is constant. Substituting from eq 143 in eq 127, we obtain 
n= 3 (145) 


A relation connecting the height at any arbitrary point with the 
partial volumes of the wave on each side of the perpendicular at this 
point is derived as follows. By virtue of the relation(1/h) (0/0zx) = 0/dz, 
when dh/dt is placed equal to zero, eq 122 yields the result: 


ae +7 ae 3) constant. 


The constant is zero, since h vanishes at infinity; therefore, 


2H? d? h ‘ 
+3 aa=0. (146) 
Integrating and employing the boundary conditions o=0, h=0; 
o=Q, h=0, we obtain the solution 


3 
h=gpp(Q—o)o 


(147) 
Since the wave is symmetrical, the maximum ordinate, /,, divides 
the wave into equal parts: n=he, c=Q/2. 
Therefore, 

3Q? 


hy ~ 16H?’ 


(148) 


which is the relation between the maximum height of the wave and 
its total volume. It was demonstrated in section IV-2 that the 
potential and kinetic energies of an intumescence are equal, to within 
a small term of the order of h*; consequently 


R=29=< (149) 


8H® 





pat TIrrotational Translation Waves 


From eq 148 and 149 we derive the following relations: 


8VE Q 8H 4 
h= 4H” and hi 39R (150) 

These relations will now be used to explain the phenomena of waves 
traveling up a channel of slowly decreasing depth. Considering eq 
85, we see that the quantity c,—c, is small, since dH/dz is small. 
Consequently, the reflected wave is negligible, and the energy of the 
transmitted wave remains sensibly constant. From eq 150 we deduce 
that hy continually increases, and Q/h; continually decreases. The 
wave gets higher and higher, and its base narrows continually. As 
a result, the wave becomes unstable and finally breaks. A quantita- 
tive analysis of this phenomenon will be given in section V—19. 


9. CNOIDAL WAVES 


In the theory developed thus far, extensive use has been made of 
the assumption that the liquid is undisturbed at infinity; that is, that 
the wave profile z=f(h) vanishes at infinity together with all its 
derivatives. In the previous section this assumption was utilized in 
deriving eq 137, and it was found that the only wave fulfilling this 
condition and propagated without change of form was the solitary 
wave. We now wish to eliminate this restriction and determine the 
form of a wave of appreciable height, governed by eq 106 and propa- 
gated without change of form. 

We therefore assume that at h=0, the derivatives 0h/Oz, 07 h/Oz’, 

. do not all necessarily vanish. The constant cz is consequently 
retained in eq 136, and we write 


by 99 LF + 3¢,=0. (151) 


Multiplying both members of eq 151 by (dn/d¢) dg and integrating: 
2 
(G2) =3bnt—39?-+3ean+ 3a, (152) 


where 3c; is another constant of integration. Since dn/d¢ is real at 
7=0, then ¢ is positive, and if 7’, 7’’, 7’’’ are the roots of the poly- 
nomial, then 


PM OPP 


nnn C1. 

Kither all the roots are positive, or one is positive and two negative, 
or one is positive and two complex. There is at least one positive value 
of n, say n’, for which dn/d¢ vanishes. If this is the only value of n for 
which the profile is horizontal, then the wave profile must become 
infinite with ¢, which is physically impossible. The third case is there- 
fore excluded. If all the roots are positive, then the wave profile must 
be included between the extreme values, since otherwise it would 
become infinite. This requires that the wave surface be everywhere 
above the primitive surface, which is likewise physically impossible, 
and the first case is also excluded. Therefore n’’ and 7’’’ are both 
negative. Writing the three roots ,,—2,—13, 732 2, we have 





82 Journal of Research of the National Bureau of Standards 


(2 2) =3(m—1) (n+) (n+n3), 


b=m—m— 1, 

Co= — ans mnt 122, 

C= 7273- 
Introducing the new variable, x, 


n=m Cos? x—1 sin? x, 


6(X)=Vvi=F sin* x, 


™ ae aa 4H 
0 seca VE” 


2. m+%M _ rth 
m+ hiths 


Selecting ¢=0 at the apex, then 77 and x=0, and we have 


eq 153 becomes 


where 








and 








(159) 


s=a[" srr make k), 


where F is an incomplete elliptic integral of the first kind. Introduc- 
ing Jacobian elliptic functions, since cos x=cn(¢/f, k), sin x= 
sn(¢/8, k) and sn? u-+cn? u=1, we have 


n=—nat (m+n) en*(§, k) (160) 


or reverting to the original coordinates 


h=—hot (uth) en'| [SA eat) fete | (161) 


This equation represents an infinite number of undulations of 
identical size and shape, each symmetrical about a vertical plane 
passing through the apex. The existence of these periodic waves 
moving without change of form was first indicated by Boussinesq 
[10, p. ~ 392], and the name, cnoidal waves, from the cn function by 
analogy with sinusoidal, is due to Korteweg and de Vries [11]. 


The relative wave length, -_ given by 
Xx al a/2 toe 
et), ae on 


F, being the complete elliptic integral of the first kind. 





me Irrotational Translation Waves 83 


The initial depth of the water in the channel being H, we must 


have 
X 


H 
ndg=0» 
or from eq 155 and eq 156 





2 (m cos? x—n2 sin? x) B 
2 : dx=0, 
J, V¥1—k sin? x . 


which from eq 158 reduces to 
(i ths) Ey (k) =hsFi(h), (163) 


where E, is the complete elliptic integral of the second kind. The 
form of the waves is completely determined by \/H, and h,/H, the 
relative wave length and height of the crest above the initial surface 
of the water, since h,/H and h;/H are determined by eq 158, 162, and 
163. 

From eq 162 we see that the wave length of the undulations be- 
comes infinite when k=1. From eq 163 if k=1, hoa=hg=0. Since 
en(x, 1)=sech x, then for k=1, eq 161 reduces to the equation of the 
solitary wave (eq 142), and the latter appears as a special case of the 
cnoidal wave. 

From the definition of } be a and from eq 154, we have 

“ A MEA 


a= +H-A- Wp (164) 


which is the expression for the velocity of the wave profile. Since 
we are now dealing with an infinite number of undulations moving in 
a channel of finite depth, a question arises as to the meaning of the 
velocity of propagation. Stokes [12, p. 202], has given two definitions 
applicable to this case, and we adopt the second one; namely, the 
velocity of propagation is the velocity of the wave profile with respect 
to the center dusts of the fluid mass included between two vertical 
planes separated by a distance equal to the wave length. 
Writing this velocity as w,: 


wo py | U+h)dz, (165) 
where U is the mean velocity in a vertical section. If U, is the mean 
velocity in the section at A=0, we have from eq 95 

U(A+h)=ho+ UoH. (166) 
Substituting in eq 165, we have for the velocity of propagation 
w@,=w—U, (167) 


== {1-2 =). (168) 


We now evaluate Uy. Integrating eq 114, we have 


oe ae OF to, 1,0, 


or 
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Up being the velocity at the bottom. It has been assumed in the 
analysis that at h=0, the surface velocity u, vanishes (see the discys- 
sion of eq 41); hence ‘from eq 114, 


H? 0 uw 
Oo Oa? 


U, _H" (52) A= 0. (169) 


From eq 104, since we can write 0/0t=—w (0/dz), we obtain by differ- 
entiating twice with respect to z, and neglecting the term containing 


O*h/ Oz": 
eae Coe =) —4 


Substituting in eq 169, we ntl have 


2 
Cenk TE. SP ge 


3 w 2 
and eq 168 becomes 
29H? /0* h 
vier 1 2 (S2),_.) 


or from eq 151, 154, 164, in terms of fy, ho, and hs, within the approxi- 
mation of the analysis, that i 3s, neglecting hyhehs/H®, 


Joel tlf —4]f1—¥].- — 


’ A=0, 


and 


40 so 60 70 80 90 


cH 
\ +0.75 


Fiaure 3.—Cnoidal waves. 
Relative displacemen of the free surface as a function of the fractional part of the relative wave length. 


This, then, is the velocity of propagation of the wave profile with 
respect to the center of inertia of the fluid mass included between 
two vertical planes separated by a distance equal to the wave length, d. 

For other treatments of this subject reference may be made to 
Boussinesq [10, p. 390], Korteweg and de Vries [11], and Lamb [4, p. 
402]. 
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Figures 3 and 4 represent two typical cnoidal waves. Figure 3 
shows on an arbitrary scale one complete undulation for two values 
of \/H, and figure 4 shows extended profiles of the same two waves 
at the correct height above the bottom, but compressed horizontally 
5times. In both cases, h,/17=0.25. 

















t=0 


Fiacure 4.—Cnoidal waves. 
Relative height above the bottom as a function of relative horizontal dimension. 


10. MOMENT OF INSTABILITY OF A SOLITARY WAVE 


Of all the intumescences of equal energy, the solitary wave is the 
only one for which the moment of instability 


@ 2 9 h3 
Mem ff) Je a2 
is & Minimum. 


The theorem may be proved by the methods of the calculus of 
variations for the solution of isometric problems. By an isometric 
problem, we mean one of the following kind: 

Among all curves joining two points (29, Ao) and (x, hy) for which 
the definite integral 


81 
K={, G(z, h, x’, h’)ds (171) 
So 
takes on a given value /, to determine the one 


Ch, r= 2(8), X= xL(S8o) 
h=K(s), ho=h(S0) 


which minimizes (or maximizes) another definite integral 


J= {. F(a, h, x’, h’)ds. 
Se 
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Here s is a parametric variable, and the primes indicate differentia. 
tion with respect to s. 
Putting 
H=F-+)4, (174) 


where \ is a numerical constant, the first necessary condition for 9 
minimum or maximum is that the curve Cy, in eq 172 must satisfy 
the differential equation 


Hy —Hiy+ Hy (2'h” —2"h’)=0, (175) 


where 
= ~—ar (176) 


The subscripts indicate partial differentiation with respect to the 
variables represented by the subscripts. The second necessary 
condition for a minimum is that 


H, 0 


at every point of the curve Cp. 

Obviously, in the present isometric problem, the integral K corre- 
sponds to the energy of an intumescence and the integral J to the 
moment of instability; that is, 


en f. "nda, (178) 
and 


=f ((8)-H- o 


We select the parametric variable s so that 

228, | (180) 
h=h(s). 
The above integrals, eq 178 and 179, become 


is { h?x'ds, 


and 


and 


Jn {Gra is, 


where a=3/H*®. Hence, 
G=h*r', 


42 
Fane oh 2’, 
z 


2 
N=s 


Mase ate 





ae an Te ae a: ae 


i tt ae ah 


Cima a: a in 
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The first necessary condition for a minimum, eq 175, after some 
simplifications, yields 


; oh 
Sah?—2Mh-+25 =0, (187) 


which also may be written as 
(188) 


if ais replaced by 3/H*. Since the left-hand member is a constant, we 
see On Comparing eq 188 with eq 112 that the intumescence having a 
maximum or & minimum moment of instability for a given energy is 
one that is propagated without change of form. The only wave which 
is propagated with constant velocity is the solitary wave. The second 
necessary condition, eq 177, is likewise satisfied; hence the moment 
of instability is a minimum, 


11. MOMENT OF INSTABILITY AND WAVE FORMATION 


The frequent occurrence of the solitary wave, owing to its produc- 
tion by a broad class of different initial disturbances, was observed 
by Scott Russell [7, p. 351; 8, p. 423]. This frequent formation of the 
solitary wave has been explained by Boussinesq [5, p. 99; 10, p. 401], 
utilizing the theorem just proved—that the solitary wave is the only 
intumescence for which the moment of instability is a minimum. An 
integration would show that this minimum is 


27 EB 


(M.)\=59 WF? (189) 


where Z is the energy of the wave and H is the depth of the channel. 

If one considers an intumescence of arbitrary form, and a solitary 
wave, both having the same energy, /, the excess of the integral M, 
of the intumescence over that of the solitary wave, M,—(M,), may 
be regarded as a measure of the rapidity with which the intumescence 
will deform during its propagation and also of the extent of the 
deformation taking place. This is why Boussinesq named the integral 
the moment of instability. 

If M,—(M,), is nner | the profile of the intumescence will oscillate 
about that of a solitary wave having the same energy, HE. The form 
of the intumescence cannot change very markedly, as this would 
require M, to increase, contrary to the theory. Friction, which has 
been ignored in the theory, will damp out these oscillations and the 
intumescence will eventually assume the shape of a solitary wave. 

It should be noted that dies does not exist a solitary wave having 
an arbitrary energy, since there is a maximum height, above which 
a solitary wave cannot be propagated without breaking (see section 
\-19). This limiting factor serves to determine the class of disturb- 
ances which can produce a solitary wave. 
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12. NEGATIVE WAVES 


Let us now examine the circumstances governing the motion of g 
negative wave limited in extent and propagated in the positive 
z-direction, negative signifying that the surface is depressed. The 
longitudinal profile of such a wave in its simplest form will consis 
at the moment of the formation of the wave of a main concave por. 
tion in the center and two convex arcs at the posterior and the anterior 
ends of the wave. As long as the form is not altered, the center of 
gravity of the wave is located at one-third of its maximum depth: 
that is, n7=h,/3. According to eq 127, the velocity of propagation of 
this wave would be 


w'=9(H—h,), (190) 


a law that is confirmed very nearly in the experiments of Bazin, 
Calculations based upon 18 observations, [6, p. 521] give as the error 
of the formula, —1.54+2.09 percent, where the + term indicates 
the standard deviation. 

The initial simple form of the wave, however, is soon modified. 
In fact, eq 112 shows that since A is negative, the third term is nega- 


Ficure 5.—Negative wave as observed by Bazin. 


tive in the concave portion and positive in the convex portion. The 
second term is more negative in the concave portion than in the con- 
vex portion, since / is smaller in absolute magnitude in the !atter 
case. Thus the head of the depression, being both at a higher eleva- 
tion than the center of the cavity and convex, while the latter is 
concave, is propagated faster than the center and becomes elongated. 
The tail of the wave also travels faster than the center for the same 
reasons and is therefore shortened as long as the requisite initial con- 
ditions obtain. However, the process of shortening cannot be indef- 
nite. Furthermore, a convex portion is necessary for the continuity 
of the surface slope; hence the convex portion is finally raised above 
the level of the primitive free surface. Since h is now positive, the 
last term of eq 112 is negative; and if either A is small enough or 
0? h/dz? large enough, it can exceed the second term, so that a condi- 
tion is reached such that w of the elevated convex portion is the same 
as w of the lowest point of the depressed body. But the elevation 
thus produced at the posterior end of the main negative body of the 
wave requires an intermediate concave portion in order to join the 
free primitive surface with a continuous slope. Such a concave por- 
tion cannot exist and cannot be formed if its velocity of propagation 
is greater than that of the positive wave preceding it. Equation 11? 
indicates that if these velocities are to be equal, the convex positive 
wave must be followed by a negative wave. Continuing the argument, 
we see that a given negative wave is followed by a series of smaller 
waves, alternatively positive and negative, and entirely convex and 
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concave, with inflection points on the level of the primitive surface 
(fig. 5). For additional figures, see [6, plate 2, fig. 3; 7, plate 52, 
fig. 10]. 


13. INTUMESCENCES OF FINITE HEIGHT BUT NEGLIGIBLE 
CURVATURE 


For an intumescence of finite height but having a free surface of 
negligible curvature, the expression of the velocity of propagation 
of volume element, eq 112, assumes the simple form 


w= VgH(1+4 7 ay (191) 


As an immediate consequence of this simplification, differential 
equation 90 can be integrated and the equation of the intumescence 
profile obtained. errs eq 191, eq 90 becomes 


ov H(1+5 Ls Hse 
which possesses the fel 


2—VgH(1 +3 bs), (193) 


where f(h) is an arbitrary function of h. The integral may be inter- 
preted as follows: If A is the elevation of the free surface at x at time 
t, then the same elevation is reproduced at z+ Az at time t+ Af; but, 
from eq 193: 


(192) 


Ar= voH(1+5 3H 


and, therefore, the intumescence height is propagated with the 
velocity 


«= VoH(1+5 7 Say (194) 


which is Airy’s expression [4, p. 243] for the velocity of propagation 
for waves of finite height. It should be noted that eq 194 represents 
the velocity of propagation of the height of an element. 

The arbitrary function, f, in eq 193 can be determined in two ways, 
depending on whether the canal is infinite in one or both directions. 
In the latter case, h is usually given for all values of x at some instant, 
t. If the relationship is in the form h=f,(z), we may invert the 
function and obtain f.(h)=2z, which is the desired function. 

In the former case, assume the channel in communication with an 
infinitely large basin and that the water level of the basin is a known 
function of time, say, A=F(t). Inverting this function, we have 
t=F\(Ah), and, if we take the point where the channel meets the basin 
as the origin of z, the arbitrary function is 


ft) =—VoHH(1+3 77 Fill). 
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In intumescences of this type, the mean velocity, U, of the fluid 
particles in a cross section is 


U=79H- H-in) (195) 


which follows from eq 96 and 191. 
14. HORIZONTAL CHANNEL IN COMMUNICATION WITH THE OCEAN 


Boussinesq has applied the method of the last article to a horizontal 
channel in communication with the ocean to determine the mean 
elevation of the ocean surface above the bed of the channel. 

We may write for the elevation of the surface water of the ocean 
above the bed of channel, 


.- Qat 
H+h=a(1+a sin i ’ z=0, (196) 


where H is the depth of water in the channel at points so remote from 
the mouth that the disturbances of the ocean are not felt, @ is the 
elevation of the mean level of the surface waters of the ocean above 
the channel bed, a the ratio of the half amplitude of the tidal oscilla- 
tions to the elevation, a, and 7 the period of the tides. Using eq 
171, we obtain for the arbitrary function, f(A), in eq 193, 


T H+h— 
I)= VoE(1+5 5 sin“ at, 


and for the free surface of water in the channel, 


2= gH 1 +3 > — sin7! Hho) (198) 


After a certain length of time, the transients will die out and the 
total volume of water entering the channel during each period 7’ will 
be zero; that is, 


(197) 


f° U(H+h)dt=0, e=0. (199) 


This condition permits the determination of H. Since U(H+h)=oh, 
eq 199 gives, substituting for w from eq 191, 


he (149 7 t=0. (200) 


Introducing h from eq 196 and integrating, we obtain 
3 (H—a)’?_ 30° 
H—a-Z Hi a Hi ’ 


which can be simplified further since H—a is a small quantity: 





H-a=3 ie. (201) 
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This is the desired relation. The mean level of the surface water 
of the ocean is below the level of the water in the horizontal channel 
at large distances from the mouth by the small amount (3/8) aa. 


15. WAVES DUE TO SUDDEN INCREASE OF DISCHARGE 


Let us examine the case in which a positive intumescence of infinite 
length is produced in a body of still water in a channel when a dis- 
charge of constant value is established at the entrance. After the 
initial disturbances have died down, the wave is of constant elevation, 
h, except in the anterior portion, and the front of the wave moves 
with a constant velocity, w. As the curvature of the layer represent- 
ing the intumescence is negligible at points sufficiently removed 
from the front, the velocity of propagation of the wave front must 


be from eq 191, 
w= o( H+5 a) (202) 
since h/H is small. 


Bazin’s experiments verify this formula closely. From his data 
(6, p. 589-541], velocities were computed and were compared with the 
theoretical values. The mean values of h and H at the two points of 
the channel between which the wave traveled when its elapsed time 
was observed were substituted in the formula. The method of 
obtaining A from the data is described on page 96. 

The data gave sufficient information to make 37 comparisons. 
Observations on eight intervals were rejected, since the observers 
reported that the wave had completely broken and reported h as zero. 

The average error of the theoretical values was —0.98+2.71 per- 
cent, where the + term indicates the standard deviation. These 
values were plotted as a function of h/H in figure 6, in an attempt to 
discover a systematic variation. The lines connect observations on 
the same wave in different intervals. The results were essentially 
negative, although there is a suggestion of an error changing in sign 
from minus to shah as h/H increases. It will be noted that the for- 
mula applies equally well to unstable waves on the point of breaking 
as well as to those having a definite h,/A ratio of 3/2. (See fig. 9.) * 

St. Venant, using the principle of momentum, has derived for the 
celerity of this type of wave, 


1 h? 
w= gH(1+5-H7+5 mp) (208) 


an expression which reduces to eq 202 when the last term is neglected. 


16. EFFECT OF CHANNEL VELOCITY DISTRIBUTION 


We may now consider a discharge wave superposed on a current 
and determine the effect of the velocity distribution in the current on 
the celerity of the wave. 

We shall assume (1) that the profile of wave does not deform dur- 
ing the motion and (2) that at the two sections, each at a distance s 
from the wave front, the velocity distributions are characteristic of 
the channel. (See fig. 7.) Let the mean velocities at these sections 
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Figure 6.—FPercentage of error of eq 202 compared with experiment, as a function 
of the relative height. 
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be U, and Us, the local velocities at any height 2 be uw, and we, and the 
depths H and H+h. We consider the body of water included at the 
instant t between the sections AB and DC. The distance between 
AB and the section containing the wave front is s; the distance be- 
tween DC and the wave front is s-+-wAt, where w is the celerity of the 
wave, and At an interval of time. The momentum of the body of 
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Figure 7.—Diagram illustrating the effect of channel velocity distribution. 
(a) Configuration at time 4; (0) Configuration at time t:+Af. 


water thus delimited may be denoted by !,+ M2, where M; is given 
by 


H 

My= post { u,dz= pwhtU,H. 
At the instant of time t+ At, the particles originally at AB have 
reached the position A’B’, and the particles originally at DC have 


reached the position D’C’. (See fig. 7.) The liquid contained 
between A’B’ and D’C’ has the momentum M,+M’,+M”, where 


H+h : 
M;=e{ (wAt—u,At)u, dz, 
0 


H 
M: =e, u’At dz. 


Introducing the Boussinesq number, n, the coefficient of the excess 
of the mean of the velocities squared over the square of the mean 
velocity, which we assume to be unaffected by depth, 
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H+h 
i: u? dz=(1+n)U? (H-+h), 


H 
f, u,? de= (1+) U2 H. 


We have, therefore, 
M,= pU,(H+h)wAt—pU/ (H+h) (1+ 0) At, 
and 
M;= p(1+n) U,? HAt. 
Thus the increase of momentum of the body of water under con- 
sideration during the time interval At is 
M,+M,'+M,’—(M,+ M2) =pdt{U,(H+h) (@— U,) + 
U,H(U,—w) —n{ U,? (H+h) — U?H}). (204) 
We may eliminate U, by the equation of continuity, 
U,(H+h) = U,H+ oh, (205) 
and obtain for the rate of change of momentum: 


ni thy 2 
a L+4!~ ol p (206 


neglecting h/H with respect to unity. 

The force which produces this rate of change of momentum is the 
difference between the total hydrostatic pressures acting on the 
sections AB and DC, which is 


rn 4s} em 


Equating 206 and 207, and assuming that the effect of the 7 term 
is small, we obtain 


(o— U)*=9H(1+5 nto)! wp } 


w= U.+\9A(1 +3 m(1-3[1- cap ]) (208) 


This is the expression for the celerity of a discharge wave or a de- 
pression wave, which takes into account the height, h, and 1, the 
channel velocity distribution. For a wave moving downstream the 

lus sign holds; for a wave moving upstream the negative sign holds. 

or & positive wave h represents an elevation and is positive; for a 
negative wave h represents a depression and is negative. 

Equation 208 can be further simplified when the approximation 
w=U,+-+79H is used. Accordingly, we obtain for a wave moving 
downstream 
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=e — 3 h n U2? 2 U, a 
mney) ie Hea) 
and for a wave moving upstream 


(210) 








— 3h [UZ 20,7 
0 Ui Vo 143 143] o-Jen) 


It is now seen that the effect of velocity distribution is appreciable 
only when the ratio U;?/gH is larger than unity, and for the ordinary 
case, in which U,?/gH is a small fraction, this effect is negligible. 


17. INITIAL WAVES 


Bazin observed that in waves of constant discharge moving in still 
water, the initial part of the wave is made up of a series of relatively 
short, alternately convex and concave, intumescences. (See fig. 8.) 


————» 
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Figure 8.—Initial wave as observed by Bazin. 


The foremost intumescence is convex, and its crest is at a higher lever 
than that of the constant elevation of the free surface of the main 
intumescence. Bazin refers to this initial part as the initial wave. 
The mechanical basis of the phenomenon was discussed by Boussinesq 
in the following manner, which is similar to the treatment in section 
V-12. 

At the head of the intumescence it is not possible for the depth to 
equal the depth of the main layer. The anterior part of the wave 
being convex, the curvature 0*h/0z? is negative and according to eq 


112, 
3h, HOh 


the velocity of propagation there is less than that part of the layer 
which follows. This layer moves with a greater velocity and thus 
overruns the anterior portion of the wave until the augmented height 
compensates for the effect of the negative curvature. This explains 
the formation of the initial wave. It is known that the height of the 
initial wave tends to oscillate around a height equal to 3/2 of the 
height, h, of the long layer. The height of the initial wave will increase 
until the velocity of propagation of its center of gravity becomes 
sensibly equal to that of the layer. But the form of this wave indi- 
cates that the center of gravity, 7, is to be placed at 1/3 of its height, 
hg; and, on the other hand, eq 127, which is applicable in the present 
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case, gives for the square of the velocity of the propagation, 9(H+,) 
Equating this to the square of the velocity of the layer, g(H +3h/2), 
eq 202, it is seen that h,=3h/2. ios 

The ratio h,/h=3/2 has been investigated by Bazin [6]. The ex. 
periments were made in a rectangular channel having a bottom slope 
of 1.5 in 1,000. The percentage change in depth being the least at 
the deep end, where the waves were generated, the motion at this end 
will nearly approximate flow in a channel of uniform depth. Caley. 
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Ficure 9.—Relative height of the head of an initial wave as a function of the relative 
depth of the main layer. 


lations made by the authors from data taken at his end of the channel 
[6, p. 542] give Ay/h=1.6140.18. Bazin [6, p. 550] considers 3/2 
to be the approximate value. 

Inasmuch as the waves are found to break as they ascend the 
channel if the original depth of the water is small enough, it is instruc- 
tive to consider this aspect of the phenomenon. Observations were 
taken on the wave at four points in the channel. The complete 
data [6, p. 539-541] have been employed to plot h,/h as a function of 
the dimensionless ratio h/H, the relative height of the wave layer. 
The average value of h was utilized in order to approximate the 
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steady height of the layer, eliminating the undulations behind the 
head of the wave. These data are represented in figure 9. The 
circles represent observations on waves which have no tendency to 
break; the squares, those in which breaking is imminent; the triangles, 
those in which breaking has occurred. An average of the 33 points 
with circles gives h,/h=1.51+0.17. Although the data do not pro- 
vide a clearly defined function, the following conclusions seem justified. 

1. As would be expected, for increasing h/H the waves become 
unstable and eventually break. The critical value is about 0.275. 

2. The value h,/h=1.5 is a reasonable approximation in the stable 
region. 

3. The ratio h,/h definitely decreases with increasing h/H, outside 
the stable region. 

Recent experiments [2, p. 158-159] indicate that h,/h is a function 
of h/H in the stable region, increasing from 1 to 2 at the critical 
value 0.28, then decreasing asymptotically toward unity in the 
unstable region. 

Where the initial wave joins the main layer, the condition of con- 
tinuity of slope requires a concave curve, that is, a curve for which 
0*h/dz? is positive. The positive curvature would entail a velocit 
of propagation greater than that of the layer as a whole, unless A 1s 
less than in the remainder of the layer. Hence a depression is pro- 
duced immediately following the initial wave. 

Continuing the argument, it is seen that the initial wave is followed 
by convexities and concavities all above the undisturbed surface. 
Because of the internal friction of the liquid, the heights of these 
convexities and concavities decrease in the posterior direction until 
they become insensible. The inflexion points where 0°h/dz? vanish 
are in a plane containing the free surface of the main layer. 


18. FAVRE’S EXPERIMENTS ON THE UNDULATIONS IN A POSITIVE 
SURGE 


In 1934, Favre [2] made an elaborate series of experiments on trans- 
lation waves. Several of these experiments were undertaken to 


arecieey the undulations a the head of a discharge wave 


[2, ch. [V]. The experiments established the following facts: 

1. The undulations are not formed immediately, but require a 
certain amount of time to establish themselves. 
= After the undulations are formed, they are of similar size and 
shape. 

3. The final configuration is a stable one. 

4. The height of the undulation at z is independent of 0, the time 
in which the flow is established, for a considerable range of values of 
z, provided Ow/z<0.2. In other words, within certain limits, a dis- 
charge which is not created instantaneously will behave like a theo- 
retically sudden increase of flow. 

5. The wave length of the undulations depends only on the mean 
height of the layer above the primitive surface. The former decreases 
as the latter increases. 

These results suggest that the undulations composing the head of 
a discharge wave ove a very general character, in the same sense 
that the solitary wave has a general character. Within certain limits, 
they are formed independently of the rate at which the discharge is 
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increased, and depend only on the final discharge rate; they assume 
a stable configuration after a short initial interval, that is, they are 
propagated without change of form. In section V-9, it was mathe. 
matically demonstrated that the only wave propagated without change 
of form is the cnoidal wave, which includes the solitary wave as q 
special limiting case. It is also true that the solitary wave is a special 
case of the discharge wave from the physical point of view. The 
solitary wave may be formed by a sudden increase of discharge which 
soon decreases to zero, while the discharge wave is caused by a sudden 
increase of discharge which is maintained. All these considerations 
inevitably suggest the question: Are the undulations composing the 
head of a discharge wave cnoidal waves? 

A comparison of figure 4 with Favre’s figure 48 |2, p. 156] supports 
an affirmative answer. The wave profiles are very similar. However, 
there is a more significant comparison which may be made, namely, 
one between the theoretical velocity of propagation, ,, and that 
observed by experiment. This is accomplished by using eq 167 and 
170. From the values of m and d/H given by the experimental data, 
we determine 3, k, and 7, by eq 163, 162, and 158, and then calculate 
w, from eq 170. Equation 191 written in the form 


w= VoH'(1 +3 7) (191) 


is used tocomputew. This equation, in the form of eq 202, has already 
been verified in section V-15. In order to calculate Up, we write 


wh’ = U,(H’ +h’), (211) 


where h’ is the average height of the layer, and H’ is the undisturbed 
depth of the liquid. It should be noted that this H’ is different from 
the H in eq 170. U, is computed from eq 211 with the assistance of 
eq 191. From eq 167 we have w,.=w—U,. These two values are 
compared in the following tables for nine nonbreaking waves. Table 1 
gives the experimental data from Favre [2, p. 157], and table 2 gives 
the results of the calculations. The explanation of the symbols is 
clearly indicated in figure 10. 
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Ficure 10.—Diagram illustrating the notation used in the studies of the positive 
surge. 
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TABLE 1.—Favre’s experimental data on positive surges 
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TaBLE 2.—Velocity comparisons between cnoidal waves and positive surges 
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The average percentage error of w, as compared to w—U, is +0.09 
+1.53, where the + term indicates the standard deviation. This 
is further evidence that the undulations in the head of a discharge 
wave are cnoidal waves. Inasmuch as these experiments were not 
intended to shed light on this particular problem, other points of 
verification cannot be easily attempted. An experiment made especi- 
ally to investigate this problem would be more suitable. 


19. THEORY OF BREAKERS 


As was indicated at the end of section V-8, there is a limiting 
height for a solitary wave, dependent on the depth of the liquid in 
the channel. If an attempt is made to form a higher wave, it breaks 
at the crest. If w is the velocity of propagation of the solitary wave, 
and u; is the absolute velocity of the liquid particle at the apex, that 
is, at the point A=/,, it is evident that the limiting height is repre- 
sented by the relation 


w—w=0, (212) 


since if u>w, the particle at the crest of the wave is traveling faster 
than the wave itself, and consequently the wave is breaking. 
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Equating eq 112 and 143 and neglecting terms in (h,/H)? or higher 
powers, we obtain for the curvature at the apex of the wave 


(3) 
Ox ? h=h 


Substituting this value in eq 116, putting z=H-+h, and neglecting 
terms of the order of (h,/H)*, the resulting velocity of the particle js 


wale XY en 


Putting this value and that from eq 143 in eq 212, the limiting 
height is a root of the equation 


hy, 3/r\ h 
At3( pH) -/ 1+ (214) 


The only positive real root of this equation is h,/H=0.731, and this 
is the limiting height of the solitary wave. 

This figure is confirmed reasonably well by experiment. Calcula- 
tions made by the authors from Bazin’s data [6, p. 510-514] give 0.71. 
Scott Russell [7, p. 352; 8, p. 425, 426, 445, 450] surmised that the 
wave broke when its height was approximately equal to the depth, 
and McCowan [13, p. 56] states that for his experiments 0.75 is a 
closer approximation for the relative elevation at the breaking point. 

The above result may now be applied to the ultimate breaking 
of waves rolling in on a gently sloping bottom. Assuming that the 
waves are solitary, we obtain from eq 148 the following limiting 
condition involving the volume Q: 


3Q? hy _ 
16H Wo": 


H=0.71yQ. (215) 


This gives the depth at which the waves will break as a function of 
the volume. If the waves travel 5 a channel with a — bed, the 


larger waves break first, and the depths at which they break is pro- 


portional to the square root of their volume. 


The authors express their appreciation to Herbert N. Eaton of the 
National Bureau of Standards for his valuable criticisms and review 
of the paper. 
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